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PREFACE 


This  study  was  conducted  at  the  Coastal  Engineering  Research  Center 
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Larson,  Assistant  Professor,  DWRE,  UL;  Dr.  Nicholas  C.  Kraus,  Senior  Scien¬ 
tist,  Research  Division  (RD) ,  CERC;  and  Dr.  Mark  R.  Byrnes,  Research  Physical 
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SBEACH:  NUMERICAL  MODEL  FOR  SIMULATING  STORM-INDUCED  BEACH  CHANGE 
NUMERICAL  FORMULATION  AND  MODEL  TESTS 

PART  I :  INTRODUCTION 

Problem  Statement  and  Objectives 

1.  Beach  stabilization  and  coastal  flood  protection  are  major  areas  of 
concern  in  the  field  of  coastal  engineering.  Beach  erosion,  accretion,  and 
change  in  the  offshore  bottom  topography  occur  naturally  through  the  transport 
of  sediment  by  waves  and  currents.  Additional  changes  result  from  perturba¬ 
tions  introduced  by  coastal  structures,  beach  fills,  and  other  engineering 
activities.  Beach  change  is  controlled  by  wind,  waves,  currents,  water  level, 
sediment  characteristics  and  supply,  and  constraints  on  sediment  movement  such 
as  those  imposed  by  coastal  structures.  These  controlling  factors  are  non¬ 
linear  and  have  great  variability  in  space  and  time.  Although  it  is  a 
challenging  problem  to  predict  the  course  of  beach  change,  such  estimations 
are  necessary  to  design  and  maintain  shore  protection  projects. 

2.  The  beach  can  be  considered  as  a  flexible  structure  that  provides 
protection  to  life  and  resources  along  the  coast.  However,  over  an  interval 
of  just  a  few  hours,  storms  can  produce  serious  damage  and  life-threatening 
situations  by  rapid  erosion  of  the  beach  and  upland  inundation.  Quantitative 
estimation  of  storm-induced  beach  change  (e.g.,  contour  retreat  and  volume  of 
eroded  material)  as  a  function  of  the  time -dependent  forcing  conditions  is 
therefore  of  great  importance.  Another  important  aspect  of  beach  response  to 
storm  action  is  the  recovery  phase  normally  occurring  at  the  end  of  a  storm 
and  during  periods  of  mild  wave  conditions.  As  wave  height  decreases  during 
the  waning  stages  of  a  storm,  waves  transform  from  erosive  to  accretive,  and 
the  foreshore  is  rebuilt.  It  is  vital  to  describe  the  process  of  recovery  if 
storm  impact  is  to  be  assessed  and  the  amount  of  material  lost  from  the  beach 
quantified.  The  recovery  process  must  also  be  modeled  if  the  impacts  of 
successive  storms  are  to  be  estimated. 

3.  Seasonal  changes  in  profile  shape  are  readily  apparent,  produced  by 
smaller,  longer  period  waves  in  the  summer  and  higher,  shorter  period  waves  in 
the  winter.  Waves  associated  with  summer  conditions  tend  to  produce  sand 
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movement  directed  toward  shore  to  build  up  the  subaerial  beach,  whereas  winter 
waves  typically  erode  the  inshore  part  of  the  beach  profile  and  deposit 
material  offshore,  creating  one  or  more  longshore  bars.  It  is  important  to 
reproduce  these  seasonal  changes  in  profile  shape  if  predictions  or  simula¬ 
tions  are  to  be  valid  for  long  periods.  Also,  the  capability  to  estimate 
variability  in  shoreline  position  due  to  seasonal  and  shorter  time  changes  in 
wave  climate  provides  a  basis  for  managing  coastal  resources  rationally  and 
identifying  if  and  when  a  net  loss  in  beach  volume  has  taken  place. 

4.  Initial  adjustments  of  beach  fill  to  wave  action  and  the  response  of 
beach  fills  to  storm  processes  are  of  central  interest  in  the  design  of  fill 
schemes  for  artificial  nourishment.  Prediction  of  protective  dune  and  beach 
profile  evolution  for  different  geometric  fill  configurations  and  grain  sizes 
can  provide  guidance  for  an  efficient  benefit-cost  analysis  of  proposed 
nourishment  designs.  Different  nourishment  templates  must  be  evaluated  and 
compared  objectively  to  assist  in  finding  the  optimal  design  to  satisfy 
criteria  based  on  storm  protection  and  recreation  benefits. 

5.  Structures  located  in  the  nearshore  zone  can  alter  beach  evolution, 
restricting  sand  transport  in  a  manner  dependent  on  the  type,  size,  and 
location  of  the  structure.  In  particular,  potential  profile  change  on  a 
seawall -fronting  beach  must  be  estimated  to  determine  if  the  wall  will  be 
damaged  during  storm  attack.  Likewise,  the  influence  of  the  wall  on  modifying 
the  beach  is  of  interest  in  project  design.  Knowledge  of  beach  change  induced 
by  seawalls  will  help  in  their  design  for  optimal  performance  of  the  structure 
and  to  appraise  its  impact  on  the  beach  (cf.  Kraus  1988). 

6.  Prediction  of  beach  evolution  with  numerical  models  has  proven  to  be 
a  powerful  technique  that  can  be  applied  to  assist  in  determining  project 
design.  Models  provide  a  framework  for  developing  project  problem  formulation 
and  solution  statements,  organizing  data  collection  and  analysis,  evaluating 
alternative  designs,  and  optimizing  the  selected  design.  Problems  related  to 
beach  profile  response  during  storms,  initial  adjustments  of  a  beach  fill, 
seasonal  changes  in  profile  shape,  and  the  influence  of  seawalls  on  profile 
evolution  are  conveniently  evaluated  using  numerical  models. 

7.  This  report  is  the  second  in  a  series  describing  development  of  a 
numerical  simulation  model  for  engineering  use  for  predicting  beach  profile 
evolution,  including  the  formation  and  movement  of  major  morphologic  features 
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such  as  longshore  bars,  troughs,  and  berms.  The  model  is  called  SBEACH,  an 
acronym  for  £torm- induced  BEAch  CHange .  The  first  report  in  the  SBEACH  series 
(Larson  and  Kraus  1989a) ,  hereafter  referred  to  as  Report  1 ,  describes  model 
development,  including  the  approach,  empirical  foundation  of  the  model,  and 
initial  tests  and  verification  with  large-scale  laboratory  data  and  limited 
field  data.  The  objectives  of  the  present  report  are  to  describe  the  numeri¬ 
cal  solution  procedure  and  to  demonstrate  applicability  of  the  model  with 
newly  acquired  field  data  pertaining  to  realistic  situations  encountered  in 
engineering  practice.  Report  3  in  this  series  will  provide  details  on  the 
operation  of  SBEACH  in  the  form  of  a  User's  Manual. 

Review  of  Methodology 

8.  The  model  SBEACH  was  developed  using  data  obtained  from  experiments 
performed  in  large  wave  tanks  with  monochromatic  waves  of  prototype -scale 
heights  and  periods.  Report  1  gives  a  thorough  discussion  of  the  data  base, 
including  a  wide-ranging  vimpirical  analysis  of  the  data  to  establish  cause - 
and-effect  relationships  between  the  waves  and  beach  profile  change.  Based  on 
this  analysis  and  physical  considerations,  a  numerical  model  was  developed  and 
tested  for  simulating  profile  evolution  generated  in  the  large  wave  tanks. 

The  model  was  then  tested  using  field  data  obtained  at  the  Coastal  Engineering 
Research  Center's  (CERC's)  Field  Research  Facility  (FRF)  in  Duck,  NC,  and  was 
shown  to  give  reasonable  agreement  between  calculated  and  measured  bar 
formation  and  movement. 

9.  In  this  report,  the  capability  of  the  model  to  predict  berm  and  dune- 
erosion  is  evaluated  using  recently  acquired  field  data  from  both  the  US  east 
and  west  coasts.  Hypothetical  cases  are  also  simulated  to  demonstrate 
applicability  and  potential  uses  of  the  model  to  predict  the  initial  adjust¬ 
ment  of  a  beach  fill  and  its  response  to  storm  action,  including  poststorm 
recovery.  Two  different  synthetic  storms  are  developed,  corresponding  to 
tropical  (hurricane)  and  extratropical  (northeaster)  events. 

10.  Improvements  in  the  numerical  model  have  been  made  beyond  those 
presented  in  Report  1.  The  wave  transformation  model  has  been  generalized  by 
using  complete  linear  wave  theory  everywhere  along  the  profile  without  any 
shallow-water  wave  approximations .  This  capability  becomes  important  in 
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applications  encountered  in  engineering  practice,  in  which  differing  wave 
conditions  can  place  the  active  calculation  at  various  locations  on  the 
calculation  grid.  Also,  it  is  now  possible  to  specify  an  arbitrary  angle 
between  the  incident  waves  and  the  bottom  contours,  thus  taking  into  account 
wave  refraction.  An  option  has  been  included  to  randomize  the  wave  input  to 
better  represent  forcing  conditions  in  the  field.  This  is  done  by  randomly 
changing  the  input  wave  height  by  a  predetermined  percentage  around  its  input 
value.  Such  a  modification  of  the  wave  height  makes  the  break  point  move  back 
and  forth  and  thus  produces  smoother  profile  features  more  in  accordance  with 
the  topography  observed  in  the  field.  Wave  runup  is  presumably  better 
represented  in  this  process,  thereby  taking  another  step  toward  improved 
representation  of  sediment  movement  on  the  foreshore. 


11.  Part  I  provides  an  introduction  and  short  review  of  Report  1. 

Part  II  reviews  the  empirical  foundation  of  the  model,  with  information  on  the 
model's  capabilities  and  limitations,  as  well  as  related  subjects  for  future 
work.  Part  III  describes  the  wave  model  that  provides  input  to  the  empirically 
developed  transport  formulas  and  a  realistic  estimate  of  the  water  level  under 
storm  surge,  tide,  and  waves,  as  described  in  Part  II.  Part  IV  describes  the 
profile  change  model,  focusing  on  the  numerical  solution  scheme  and  data  input 
requirements.  Part  V  first  reviews  sensitivity  tests  of  the  model  and  then 
describes  two  case  studies,  one  for  the  US  east  coast  and  one  for  the  US  west 
coast.  In  Part  VI,  the  model  is  applied  for  hypothetical  examples  to  demon¬ 
strate  a  range  of  uses  in  beach-fill  design.  Part  VII  provides  concluding 
observations  of  the  present  capabilities  and  limitations  of  SBEACH,  leading  to 
recommendations  for  further  study  and  development  of  the  model. 
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PART  II:  REVIEW  OF  EMPIRICAL  FOUNDATION 


12.  This  chapter  gives  a  summary  of  the  empirical  foundation  of  the 
beach  profile  change  numerical  simulation  model  SBEACH.  The  data  base  used  in 
model  development  is  outlined,  and  the  relevance  of  laboratory  data  is  dis¬ 
cussed  with  reference  to  field  conditions.  Salient  aspects  of  cross -shore 
sand  transport  characteristics  derived  from  the  data  are  presented,  focusing 
on  direction  and  distribution  of  transport.  The  four  different  transport 
regions  used  in  the  model  are  introduced  and  related  to  findings  in  nearshore 
wave  dynamics.  In  addition,  the  semi -empirical  cross-shore  sand  transport 
rate  equations  are  summarized. 


Data  Base 


13.  The  complexity  and  randomness  of  nearshore  sediment  processes,  the 
limited  capability  of  quantifying  nearshore  fluid  and  sediment  movement,  and 
the  limited  availability  of  input  data  in  realistic  applications  all  indicate 
that  numerical  models  for  simulating  beach  profile  change  in  engineering 
situations  must  be  based  on  empirical  relations.  In  Report  1,  it  was  argued 
that  a  macroscale  approach  of  model  development  based  on  a  firm  empirical 
foundation  provided  a  reasonable  approach  which  paralleled  the  regularities  in 
morphologic  change  exhibited  by  beaches.  These  changes  are  on  a  spatial  scale 
of  meters  and  on  a  time  scale  of  hours  to  days.  Such  regularity  contrasts 
markedly  from  the  randomness  and  complexity  of  individual  particle  movement. 
However,  the  lack  of  detailed  field  data  encompassing  high- resolution  measure¬ 
ments  in  time  and  space  together  with  associated  wave  and  water  level  condi¬ 
tions  made  it  necessary  to  use  laboratory  data  to  enhance  development  of  such 
a  model. 

14.  Use  of  laboratory  data  must  address  the  problem  of  scaling  in 
attempts  to  transfer  observations  and  results  to  the  prototype.  It  is  desir¬ 
able  to  perform  experiments  at  a  scale  which  is  close  to  that  of  waves 
producing  significant  change  on  a  natural  beach,  without  distorting  the  kine¬ 
matic  and  dynamic  interaction  between  fluid  and  sediment  motion.  However, 
only  a  small  number  of  experiment  programs  have  been  carried  out  on  a  scale 
comparable  to  the  field  situation  (Saville  1957,  Kajima  et  al.  1982,  Vellinga 
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1982,  Dette  and  Uliczka  1987).  Data  from  the  first  two  laboratory  programs 
were  used  for  numerical  model  development  in  this  study  (Report  1)  because 
they  provided  an  extensive  range  of  conditions  that  reasonably  replicated  a 
naturally  sloping  beach.  In  the  combined  data  set,  a  large  number  of  wave  and 
beach  conditions  were  available  for  analysis,  and  the  empirical  relations 
obtained  were  less  likely  to  be  influenced  by  experiment-specific  conditions. 

15.  Model  development  was  largely  based  on  data  obtained  in  the 
prototype -scale  experiments,  whereas  model  verification  also  used  data  from  a 
natural  beach.  Larson  (1988)  and  Larson  and  Kraus  (1989a)  give  a  detailed 
discussion  of  the  two  laboratory  experiment  programs  together  with  the  results 
of  extensive  data  analysis  of  the  geometric  characteristics  of  the  profile  and 
their  relation  to  wave  and  sand-beach  properties. 

16.  In  this  report,  only  a  short  review  of  the  experiments  is  presented 
to  provide  a  background  for  discussion  of  the  numerical  model.  The  first  data 
set  was  obtained  in  experiments  performed  by  the  US  Army  Corps  of  Engineers 
(CE)  in  the  years  1956-57  and  1962  (Saville  1957,  Caldwell  1959,  Kraus  and 
Larson  1988)  at  Dalecarlia  Reservation,  Washington,  DC.  The  second  data  set 
was  derived  from  experiments  carried  out  at  the  Central  Research  Institute  of 
Electric  Power  Industry  (CRIEPI)  in  Japan  (Kajima  et  al.  1982,  1983). 

CE  experiments 

17.  The  CE  experiments  were  conducted  in  a  concrete  tank  193.5  m  long, 
4.6  m  wide,  and  6.1  m  deep  and  having  a  standard  operating  depth  of  4.6  m  (for 
details  of  the  experimental  facility  see  Raynor  and  Simmons  (1964)  and  Kraus 
and  Larson  (1988)).  Most  cases  were  started  from  a  plane  initial  beach  slope 
of  1:15,  and  waves  of  constant  height  and  period  were  applied  until  near¬ 
equilibrium  conditions  were  attained,  A  run  typically  lasted  40  to  60  hr,  and 
the  water  level  was  held  constant  for  all  cases  except  one  in  which  a  step- 
type  sinusoidal  variation  was  employed.  On  the  order  of  10  to  15  profile 
surveys  were  made  during  a  typical  case,  and  the  survey  interval  was  1.2  m. 

18.  Two  different  sands  were  used  with  median  grain  sizes  of  0.22  and 
0.40  mm.  Eighteen  distinct  cases  have  been  documented  (Kraus  and  Larson  1988) 
of  which  sixteen  started  from  a  plane  slope.  Wave  height  varied  between  0.55 
and  1.68  m  in  the  horizontal  portion  of  the  tank,  and  wave  period  ranged 
between  3.75  and  16.0  sec.  Some  observations  were  made  of  the  breaking  wave 
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height  and  location  during  the  runs,  and  a  limited  number  of  wave  runup 
measurements  were  performed  as  well. 

CRIEPI  experiments 

19.  The  CRIEPI  experiments  were  similar  to  the  CE  experiments  in 
design,  and  wave  tank  dimensions  (205  by  3.4  by  6  m)  were  comparable  to  those 
of  the  CE  tank.  The  experiment  program  consisted  of  24  cases  with  a  varying 
initial  beach  profile,  including  both  plane  initial  slopes  and  profiles  of 
arbitrary  shape.  The  generated  wave  height  ranged  between  0.3  and  1.8  m  and 
the  wave  period  between  3  and  12  sec.  All  cases  were  run  with  monochromatic 
waves  and  a  fixed  water  level.  For  cases  starting  from  a  plane  beach  slope 
(17  of  the  24  cases),  the  initial  slope  varied  from  1:50  to  1:10.  Two 
different  sands  with  median  diameters  of  0.27  and  0.47  mm  were  used. 

20.  Wave  height  and  wave  setup  along  the  profile  were  measured.  The 
number  of  profile  surveys  made  during  a  run  was  somewhat  fewer  than  for  the  CE 
experiments,  but  the  spatial  resolution  of  the  measurements  was  considerably 
higher.  Refer  to  Report  1  for  details  of  data  analysis. 

Relevance  of  large  wave  tank  experiments 

21.  Although  scale  effects  in  the  large  wave  tanks  were  probably  small, 
some  differences  in  response  of  the  beach  profile  as  compared  with  a  natural 
beach  are  expected.  The  main  questions  concern  the  relevance  of  the  labora¬ 
tory  data;  that  is,  have  the  main  mechanisms  governing  cross -shore  sand 
transport  on  a  natural  beach  been  reproduced  in  the  wave  tank,  and  have 
additional  effects  been  induced  in  the  tank  environment  which  are  not  found  on 
a  beach?  One  procedure  for  evaluating  the  appropriateness  of  wave  tank  data 
is  to  compare  observed  profile  evolution  with  that  found  in  the  field;  such 
comparisons  are  summarized  below. 

22.  The  speed  of  bar  movement  calculated  from  the  large  wave  tank  data 
(Report  1)  was  very  similar  to  observations  made  by  Sallenger,  Holman,  and 
Birkemeier  (1985)  on  a  natural  beach  during  a  storm  event.  In  general,  the 
time  scale  of  profile  change,  both  onshore  and  offshore,  was  rapid  for  the 
field  measurements  and  in  accordance  with  what  was  observed  in  the  large  wave 
tank  experiments.  In  particular,  onshore  bar  movement  occurring  during  the 
poststorm  recovery  stage  took  place  rapidly  in  the  field  just  as  was  noted 
during  the  accretionary  cases  in  the  large  wave  tanks. 
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23.  The  ratio  between  depth  to  trough  bottom  and  depth  to  bar  crest 
obtained  in  the  large  tanks  displayed  properties  in  good  agreement  with  field 
observations  reported  by  Shepard  (1950).  Bars  in  the  field  exposed  to  waves 
with  longer  wave  periods  showed  a  smaller  value  of  the  ratio  between  depth  to 
trough  and  depth  to  crest  in  accordance  with  the  wave  tank  data.  The  value  of 
this  ratio  was  preserved  during  offshore  bar  movement  in  the  large  wave  tank 
experiments;  Sallenger,  Holman,  and  Birkemeier  (1985)  also  found  this  to  be 
the  case  during  a  storm  event  recorded  in  the  field. 

24.  A  criterion  to  predict  overall  profile  response  (on-  or  offshore 
sand  movement)  based  on  the  sand  fall  speed  and  the  deepwater  wave  steepness 
was  derived  from  the  large  wave  tank  data.  Tnis  criterion  proved  to  be  valid 
for  natural  beaches  if  the  mean  wave  height  was  used  in  the  criterion  as  the 
representative  statistical  wave  height.  The  two  empirical  coefficients  in  the 
formula  expressing  the  criterion  were  the  same  for  the  wave  tank  data  and 
field  data,  further  emphasizing  the  relevance  of  wave  tank  experiments  as  well 
as  the  apparent  lack  of  scale  effects  in  the  experiments. 

25.  Hallermeier  (1987)  pointed  out  the  pertinence  of  results  obtained 
in  the  large-scale  laboratory  tests.  He  found  quantitative  agreement  in  terms 
of  shoreline  retreat  and  erosive  geometry  between  a  large  wave  tank  case  and 
measured  erosion  during  a  storm  at  Bethany  Beach,  Delaware,  where  the  condi¬ 
tions  in  the  tank  and  in  the  field  were  similar.  The  capability  of  SBEACH  to 
predict  eroded  volume  as  well  as  breakpoint  bar  movement  in  field  situations 
also  supports  the  validity  of  large  wave  tank  studies  (Report  1) . 

26.  An  obvious  restriction  of  the  wave  tank  data,  which  is  an  advantage 
from  the  viewpoint  of  clarifying  cross -shore  transport  characteristics,  is  the 
absence  of  longshore  fluid  and  sand  motion.  However,  one  of  the  main  reasons 
for  conducting  profile  change  experiments  in  wave  tanks  is  to  gain  insight 
into  the  physical  processes  of  cross -shore  sand  movement.  The  large  wave  tank 
experiments  were  conducted  using  waves  with  a  constant  height  and  constant 
period  (monochromatic  waves)  during  a  run.  Thus,  it  is  an  underlying  assump¬ 
tion  that  quantitative  understanding  of  profile  evolution  may  be  obtained  by 
studying  the  response  to  monochromatic  waves  of  uniform  height.  In  the  model, 
the  effect  of  irregular  waves  is  simulated  by  using  a  time  series  of  waves  of 
different  heights  and  periods.  This  procedure  is  probably  adequate  for 
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simulating  profile  change  by  individual  irregular  waves  but  neglects  the 
contribution  of  other  types  of  waves  such  as  infragravity  waves. 

27.  An  alternative  approach  is  to  employ  a  statistical  or  probabilistic 
description  of  the  incident  waves,  which  would  certainly  lead  to,  for  example, 
more  gently  sloping  bars  that  develop  in  the  field  as  opposed  to  the  steeper 
bars  that  develop  under  monochromatic  waves  of  constant  height  in  tanks.  A 
proper  limit  of  such  models  is  replication  of  profile  change  and  bar  growth 
under  an  extremely  narrow-banded  wave  height  and  period  distribution.  In 
development  of  SBEACH,  there  was  concern  that  an  a  priori  probabilistic  or 
statistical  description  would  introduce  smoothing  and  mask  or  obliterate  the 
significance  of  individual  wave  components.  Therefore,  the  approach  used  is 
founded  on  profile  change  produced  under  elemental  monochromatic  waves  of 
constant  height  and  proceeds  to  describe  change  under  irregular  waves  by 
superposition. 

28.  Randomly  varying  wave  height,  period,  and  direction  produce  a  break 
point  that  continually  changes  position,  leading  to  smoother  profile  develop¬ 
ment  with  features  less  pronounced  than  in  the  monochromatic  wave  case 
(Keulegan  1948).  In  the  large  wave  tank  experiments,  beach  slopes  tended  to 
become  steep  as  the  profile  approached  equilibrium,  especially  around  the  bar 
or  at  the  foreshore  step  (Report  1) .  This  was  attributed  to  the  applied 
constant  wave  height,  wave  period,  and  water  level,  which  promoted  gradual 
breakpoint  movement  and  the  consequent  buildup  of  steep  beach  slopes.  In 
summary,  a  fundamental  assumption  guiding  model  development  is  that  short- 
period  breaking  waves  during  storm  events  are  the  dominant  cause  of  beach 
change  in  the  surf  zone  and  thus  of  bar  formation  and  movement. 

Long-period  waves 

29.  Investigations  on  several  coasts  (e.g.,  Guza  and  Thornton  1982, 
Greenwood  and  Sherman  1984,  Wright  and  Short  1984)  have  found  that  significant 
amounts  of  energy  can  be  contained  in  lower  frequencies  of  the  wave  power 
spectra  in  very  shallow  water  and  in  the  swash.  The  presence  of  lower 
frequency  wave  motion  cannot  be  predicted  at  present,  but  it  appears  to  be 
significant  in  association  with  storms.  It  has  also  been  shown  that  long 
waves  and  periodic  arrival  of  high  (short-period)  wave  groups  can  be  respon¬ 
sible  for  as  much  or  more  sediment  suspension  than  individual  short-period 
waves  (e.g.,  Downing  1984;  Sternberg,  Shi,  and  Downing  1984;  Beach  and 
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Sternberg  1987)  in  shallow  portions  of  the  surf  zone.  Thus,  if  long-period 
waves  or  wave  groups  are  present,  their  contribution  to  the  total  sediment 
transport  is  expected  to  be  significant.  An  equally  if  not  more  important 
role  played  by  long  waves  and  wave  groups  than  in  direct  movement  of  sediment 
is  their  contribution  to  elevating  water  level  near  the  shore.  This  allows 
short-period  waves  to  attack  higher  on  the  beach  than  might  be  expected. 

30.  Quantitative  engineering  relationships  for  predicting  sediment 
transport  produced  by  long  waves  and  groups  have  not  yet  been  developed. 
Therefore,  their  contribution  is  not  directly  represented  in  the  model.  When 
such  transport  relations  become  available,  it  should  be  possible  to  incor¬ 
porate  them  in  SBEACH  if  they  can  be  expressed  in  a  macroscale  framework. 

Field  data 

31.  Field  data  from  the  FRF  ware  used  to  verify  the  numerical  model, 
primarily  with  regard  to  bar  movement  (Report  1).  Temporal  changes  in 
shoreline  position  at  the  FRF  are  relatively  small,  making  the  data  poten¬ 
tially  unsuitable  for  testing  the  predictive  capabilities  of  the  model  for 
evaluating  dune  and  foreshore  change.  In  the  present  report,  dune  and 
foreshore  change  predictions  are  evaluated  using  two  recently  acquired  field 
data  sets  augmenting  field  verification  given  in  Report  1. 

Cross -Shore  Transport  Direction  and  Distribution 

32.  Cross -shore  movement  of  sediment  is  governed  by  properties  of  the 
velocity  field  and  sediment  concentration.  Thus,  a  first-principles  descrip¬ 
tion  of  the  waves,  velocity,  and  concentration  is  necessary  not  only  across 
the  surf  zone  but  also  through  the  water  column  at  every  point  on  the  calcula¬ 
tion  grid.  Presently,  such  a  detailed  prediction  is  not  possible  in  numerical 
models  of  beach  profile  change  for  engineering  use  to  describe  the  evolution 
of  natural  beaches.  However,  Dally  and  Dean  (1984)  and,  more  recently, 
Roelvink  and  Stive  (1989)  have  had  success  with  first-principles  type  models 
in  reproducing  specific  occurrences  of  profile  change. 

33.  Sediment  concentration  in  the  surf  zone  is  closely  related  to  the 
generation  of  turbulent  motion,  which  depends  on  wave  breaking.  Accordingly, 
it  is  plausible  that  sediment  concentration,  and  thus  the  amount  of  material 
available  for  transport,  is  closely  related  to  wave  energy  diqsipatio 
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produced  by  wave  breaking  (Dean  1977,  Kriebel  1982,  Moore  1982,  Kriebel  and 
Dean  1985).  This  macroscale  approach  is  used  in  SBEACH.  In  the  model, 
prediction  of  the  cross-shore  velocity  field  is  avoided  by  expressing  the 
cross-shore  sand  transport  rate  directly  in  terms  of  wave  and  profile  proper¬ 
ties.  The  direction  of  transport  is  predicted  from  an  empirical  criterion 
derived  from  the  large  wave  tank  data  and  verified  with  field  data,  and  the 
magnitude  of  transport  in  the  main  portion  of  the  surf  zone  is  related  to  the 
wave  energy  dissipation  per  unit  water  volume.  This  simplified  approach 
implies  a  transport  rate  distribution  directed  either  onshore  or  offshore 
along  the  entire  profile  at  a  specific  instant  in  time. 

Cross -shore  transport  direction 

34.  Criteria  for  predicting  whether  a  beach  will  erode  or  accrete 
through  cross -shore  sand  transport  processes  have  been  suggested  by  a  number  of 
authors  (see  Larson  and  Kraus  1989a) .  Most  criteria  refer  to  general  types  of 
profile  configuration,  such  as  winter/summer,  storm/normal,  bar/berm,  dissipa¬ 
tive/reflective,  which  evolve  under  certain  wave  and  sediment  conditions, 
whereas  some  authors  have  focused  on  predicting  the  transport  direction.  In 
general,  there  is  no  contradiction  between  the  two  types  of  criteria,  and  they 
can  be  regarded  as  equivalent  in  a  macroscopic  sense.  A  profile  on  which  bars 
develop  has  transport  directed  predominantly  offshore,  whereas  a  profile 
exhibiting  berm  buildup  mainly  experiences  onshore  sand  movement.  However,  if 
a  criterion  for  predicting  transport  direction  is  expressed  in  terms  of  local 
wave  properties,  there  could  be  a  significant  difference  between  transport 
direction  and  the  overall  profile  type  that  evolves. 

35.  Criteria  for  predicting  profile  type  or  overall  transport  direction 
typically  include  one  parameter  for  characterizing  the  incident  wave  condition 
and  another  parameter  involving  some  property  of  the  sediment  (grain  size  or 
fall  speed)  and/or  beach  slope.  Deepwater  wave  steepness  HQ/L0  is  the  most 
commonly  used  descriptor  of  the  waves  (Waters  1939,  Rector  1954,  Iwagaki  and 
Noda  1962,  Nayak  1970,  Dean  1973,  Sunamura  and  Horikawa  1974),  whereas  several 
different  parameters  for  describing  sediment  properties  have  been  proposed, 
such  as  DS0/Lo  (Rector  1954),  Ho/D50  (Iwagaki  and  Noda  1962),  and  nw/gT 
(Dean  1973) ,  in  which 

For  convenience,  symbols  and  abbreviations  are  listed  in  Appendix  A. 
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HQ  -  deepwater  wave  height,  m 
LQ  -  deepwater  wavelength,  m 
D50  -  median  grain  size,  mm 
w  -  sand  fall  speed,  m/sec 
g  -  acceleration  due  to  gravity,  m/sec2 
T  -  wave  period,  sec 

In  development  of  SBEACH,  the  deepwater  wave  steepness  H0/L0  and  the  para¬ 
meter  H0/wT  ,  called  the  dimensionless  fall  speed  or  fall  speed  parameter, 
were  found  to  give  a  good  distinction  between  profiles  exhibiting  mainly  bar 
and  berm  formation.  The  fall  speed  parameter  also  emerged  naturally  in  the 
course  of  regression  analysJs  relating  morphodynamic  properties  of  the  profile 
to  the  incident  waves  (Larson  1988;  Larson,  Kraus,  and  Sunamura  1988;  Larson 
and  Kraus  1989a) . 

36.  Figure  1  shows  the  criterion  for  distinguishing  erosion  and 
accretion  developed  for  data  from  the  large  wave  tank  experiments.  Separation 
of  eroded  and  accreted  (bar  and  berm)  profiles  is  achieved  by  the  equation 
expressing  the  solid  line  drawn  in  the  figure  and  given  by 


% 

L0 
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in  which  M  -  0.00070  is  an  empirically  determined  coefficient.  If  the  left 
side  of  Equation  1  is  less  (greater)  than  the  right  side,  the  profile  is 
predicted  to  erode  (accrete).  Wave  steepness  describes  the  wave  asymmetry, 
whereas  wave  height  and  period  appearing  in  the  fall  speed  parameter  account 
for  the  absolute  magnitudes  of  those  quantities.  The  fall  speed  accounts  for 
the  settling  characteristics  of  sand  particles.  Saville  (1957)  showed  on  the 
basis  of  the  early  CE  experiments  that  the  magnitude  of  the  wave  height 
controls  erosion  and  accretion  in  addition  to  the  wave  steepness. 

37.  Viewing  Figure  1,  it  is  seen  that  higher  waves  and  shorter  periods 
tend  to  promote  erosion.  Coarser  grained  beaches,  for  which  the  fall  speed  is 
greater,  tend  to  more  readily  accrete  than  do  finer  grained  beaches  exposed  to 
the  same  wave  conditions.  These  predictions  conform  with  field  observations. 

38.  In  describing  profile  change  in  the  wave  tank  studies,  specifica¬ 
tion  of  wave  and  sediment  properties  to  be  used  in  Equation  1  is  unambiguous 


17 


1.0000 


0.1000 


0.0100 


0.0010 


0.0001 

0.1  1.0  ty/WT  10.0 

Figure  1.  Criterion  for  predicting  beach  erosion  and  accretion 
based  on  monochromatic  laboratory  wave  data 

because  of  the  fixed  wave  height  and  period  and  well-defined  sand  size. 
However,  in  order  to  apply  the  equation  to  a  natural  beach,  appropriate  wave 
and  sediment  characteristics  must  be  available.  Field  data  were  collected 
from  various  investigations  around  the  world  to  establish  the  proper  statis¬ 
tical  descriptors  to  be  used  in  Equation  1.  The  data  used  are  summarized  in 
Sunamura  (1980a)  and  Seymour  (1986),  where  movement  of  the  deeper  contours  was 
used  to  evaluate  profile  change  and  transport  direction  in  the  latter  investi¬ 
gation.  Report  1  gives  a  more  complete  description  of  the  field  data  used  in 
evaluating  the  performance  of  Equation  1. 

39.  Three  deepwater  statistical  wave  heights  (mean  wave  height  HQ  , 
root -mean- square  wave  height  Hrms  ,  significant  wave  height  Hso)  were  tested 
in  Equation  1  to  obtain  the  best  division  between  accretionary  and  erosional 
profiles  for  field  measurements.  Although  reasonable  delineation  could  be 
obtained  with  each  of  the  wave  heights,  use  of  mean  wave  height  provided  the 
best  delineation  with  the  same  value  of  M  as  found  for  the  large  wave  tank 
data.  Figure  2  displays  in  analogy  to  Figure  1  the  criterion  (solid  line) 
together  with  the  data  from  natural  beaches.  If  deepwater  significant  wave 
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Figure  2.  Criterion  for  predicting  beach  eri.  :ion  and  accretion 


height  was  reported  in  the  field  data,  mean  wave  height  was  calculated  as 
H0  -  0.626  Hso  assuming  a  Rayleigh  distribution.  The  mean  or  peak  period  of 
the  spectra  was  chosen  as  the  representative  period  to  be  used  in  Equation  1 
together  with  median  grain  size  of  the  beach  sediment. 

40.  Mimura,  Otsuka,  and  Watanabe  (1986)  compared  profile  change,  trans¬ 
port  direction,  and  transport  rate  produced  in  a  small  wave  tank  in  comparable 
cases  using  regular  and  irregular  waves.  They  observed  that  the  criterion 
developed  by  Sunamura  and  Horikawa  (1974)  gave  the  best  delineation  between 
erosion  and  accretion  if  mean  wave  height  was  used  in  analyzing  the  irregular 
wave  tests.  Profile  change  proceeded  at  a  slower  rate  for  irregular  waves, 
and  this  was  attributed  to  the  presence  of  both  accretionary  and  erosional 
wave  components  in  a  wave  train.  Their  observation  supports  the  method  of 
superimposing  the  effects  of  individual  waves  to  obtain  profile  evolution  for 
irregular  waves  in  SBEACH.  Kraus  and  Horikawa  (1990)  further  discuss  the 
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Cross-shore  sand  transport  rate  distribution 

41.  Cross -shore  sand  transport  rate  distributions  were  calculated  from 
consecutive  profiles  in  time  by  integrating  the  mass  conservation  equation. 
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This  technique  has  been  applied  by  van  Hijum  (1974,  1976),  Hattori  and 
Kawamata  (1980),  and  Shimizu  et  al.  (1985).  The  transport  rate  obtained 
constitutes  an  average  net  rate  for  the  period  between  the  two  profile 
surveys.  Report  1  describes  a  wide  range  of  analyses  of  cross -shore  transport 
rate  characteristics  and  transport  rate  formulas  derived  on  the  basis  of 
physical  considerations  and  observations  from  the  data.  In  the  following, 
some  important  aspects  of  the  cross -shore  sand  transport  rate  are  discussed 
regarding  typical  distribution  shapes  in  time  and  space. 

42.  An  important  property  of  the  cross -shore  sand  transport  rate  as 
calculated  from  measured  profile  change  is  the  temporal  decay  of  the  rate  of 
change  for  profiles  exposed  to  constant  waves  and  water  level.  The  calculated 
maximum  sand  transport  rate  was  more  than  one  order  of  magnitude  smaller  at 
the  end  of  a  run  as  compared  with  the  earliest  part  of  the  run.  A  small 
amount  of  transport  activity  still  occurred  at  the  end  of  the  runs,  attributed 
in  part  to  profile  change  introduced  by  small  variations  in  forcing  conditions 
and  randomness  in  fluid  and  sand  particle  motion.  Thus,  in  fitting  empirical 
functions  by  least-squares  to  calculated  maximum  transport  rates,  an  inverse 
dependence  on  time  gave  a  better  description  than  an  exponential  decay  because 
of  the  finite  amount  of  transport  prevailing  at  the  end  of  the  run. 

43.  By  calculating  the  cross-shore  transport  rate  distribution  using 
initial  and  final  profiles  of  a  run,  termed  equilibrium  distributions,  a 
picture  of  the  overall  profile  response  is  obtained.  Equilibrium  transport 
rate  distributions  were  classified  as  three  main  types  (Larson  and  Kraus 
1989a).  Figure  3  illustrates  the  characteristics  of  these  types  by  displaying 
the  calculated  distributions  from  selected  large  wave  tank  cases.  The  equili¬ 
brium  distribution  types  may  be  summarized  as: 

a.  Erosional  (Type  E) ;  transport  directed  offshore  along  the 
entire  profile. 

b.  Accretionary  (Type  A) ;  transport  directed  onshore  along  the 
entire  profile. 

c.  Mixed  accretionary  and  erosional  (Type  AE) ;  transport  directed 
onshore  along  the  seaward  part  of  the  profile  and  transport 
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44.  A  majority  of  cases  fell  in  the  classifications  of  Type  E  or 
Type  A,  indicating  predominantly  unidirectional  transport  during  the  run. 

Cases  having  Type  AE  distributions  correspond  to  profiles  where  little  change 
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Figure  3.  Net  cross -shore  sand  transport 
rate  distributions 
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place,  implying  a  transport  direction  that  was  not  so  distinct.  If  signi¬ 
ficant  profile  change  occurred,  transport  was  directed  primarily  onshore  or 
offshore.  Thus,  in  the  numerical  model,  which  has  as  a  primary  objective 
simulation  of  events  producing  significant  profile  change,  it  is  reasonable  to 
assume  a  unidirectional  transport  rate  distribution  as  a  first  approximation. 

45.  Sand  transport  rate  distributions  calculated  from  consecutive 
profile  surveys  did  not  show  as  clear  a  trend  as  did  the  equilibrium  distribu¬ 
tions,  particularly  for  consecutive  surveys  taken  at  later  times.  In  the 
beginning  of  a  run,  when  large  profile  change  took  place,  the  transport  rate 
distribution  was  typically  unidirectional.  Complexity  in  transport  rate 
distributions  at  later  times  probably  reflects  a  close  balance  among  the  many 
factors  contributing  to  sand  movement,  for  which  small  instabilities  away  from 
equilibrium  could  produce  relatively  large  localized  contributions  to  the 
transport.  In  order  to  model  beach  change  in  situations  where  the  beach  is  in 
near-equilibrium  with  the  water  level  and  incident  waves,  prediction  of  more 
complex  distributions  will  probably  have  to  be  included. 

Transport  Regions  and  Transport  Rates 

46.  Report  1  describes  analyses  of  individual  net  sand  transport  rate 
distributions  with  the  aim  of  relating  the  local  transport  rate  to  wave  and 
sediment  characteristics.  Research  in  nearshore  wave  dynamics  has  identified 
regions  of  different  hydrodynamic  properties  associated  with  wave  shoaling, 
wave  decay,  and  intensity  of  small-  and  large-scale  motion.  Svendsen,  Madsen, 
and  Buhr  Hansen  (1978)  defined  an  outer  and  inner  region  of  the  surf  zone 
according  to  the  character  of  wave  and  wave -induced  water  motion.  The  outer 
region,  located  immediately  shoreward  of  the  wave  break  point,  is  distin¬ 
guished  by  a  rapid  transition  in  wave  height  with  large-scale  flow  patterns. 

At  some  location  landward  of  the  break  point,  the  large-scale  flow  decays  and 
turns  into  small-scale  turbulent  fluctuations;  simultaneously  the  wave  changes 
character  and  becomes  borelike.  The  region  shoreward  of  this  point  is  known 
as  the  inner  region  and  extends  to  the  location  where  wave  runup  begins. 

47.  The  concept  of  two  different  regions  in  the  surf  zone  had  already 
been  suggested  by  Miller  (1976),  who  distinguished  a  vortex  and  a  bore  region 
(see  also  Basco  1985  and  Jansen  1986).  Miller  also  noted  that  the  main 
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difference  between  plunging  and  spilling  breakers  was  the  generation  of 
vortices  in  the  outer  region.  In  the  inner  (bore)  region,  flow  conditions 
were  similar  and  independent  of  breaker  type,  as  also  noted  by  Svendsen, 
Madsen,  and  Buhr  Hansen  (1978).  Svendsen  (1987)  analyzed  flow  measurements  to 
determine  the  turbulent  kinetic  energy  in  the  surf  zone.  He  found  that  the 
variation  over  depth  of  che  time -mean  turbulent  kinetic  energy  was  small 
(Peregrine  and  Svendsen  1978) ,  and  he  attributed  this  to  a  strong  mixing 
created  by  large-scale  vortices.  Skjelbreia  (1987)  studied  breaking  of 
solitary  waves  in  the  laboratory  and  defined  four  different  zones  of  wave 
transformation:  gradual  and  rapid  shoaling,  and  gradual  and  rapid  decay. 

48.  The  identification  of  regions  with  different  wave  characteristics 
in  the  nearshore  stimulated  a  separate  analysis  of  the  calculated  net  trans¬ 
port  rates  in  these  regions.  An  important  assumption  underlying  the  numerical 
model  SBEACH  is  the  significance  of  breaking  and  broken  waves  in  determining 
the  cross-shore  sand  transport  rate.  Quite  different  transport  character¬ 
istics  should  be  found  inside  and  outside  the  surf  zone,  making  it  reasonable 
to  assume  that  the  transport  relationships  are  different  in  a  macroscopic 
sense.  Furthermore,  the  dynamics  of  transport  in  the  swash  are  governed  by 
properties  of  the  uprush  bore  and  the  backwash,  sediment  characteristics 
(grain  size,  porosity,  grain  density,  etc.),  and  beach  slope,  which  makes  it 
logical  to  define  the  swash  zone  as  another  distinct  region  of  transport. 

49.  Based  on  the  division  of  the  profile  applied  in  nearshore  wave 
dynamics  and  the  physical  characteristics  of  sediment  transport  under  various 
flow  conditions,  four  different  zones  of  transport  (Figure  4)  were  introduced 
(Larson,  Kraus,  and  Sunamura  1988;  Larson  and  Kraus  1989a).  These  zones  are: 

a.  Zone  I:  From  the  seaward  depth  of  effective  sand  transport 

to  the  break  point  (prebreaking  zone) . 

b.  Zone  II:  From  the  break  point  to  the  plunge  point  (breaker 

transition  zone) . 

c.  Zone  III:  From  the  plunge  point  to  the  point  of  wave  reforma¬ 
tion  or  to  the  swash  zone  (broken  wave  zone) . 

d.  Zone  IV:  From  the  shoreward  boundary  of  the  surf  zone  to  the 
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50.  The  zone  located  between  the  break  point  and  the  plunge  point  is 
analogous  to  the  outer  region  (vortex  region)  proposed  by  Svendsen,  Madsen, 
and  Buhr  Hansen  (1978).  A  certain  distance  is  required  after  breaking  before 
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where 


Zone  IV: 


‘x  -  xr1 

q  ■  q*£-ny 


Xr  <  X  <  Xz 


(5) 


q  -  net  cross -shore  sand  transport  rate,  m3/m-sec 
A1(2  -  spatial  decay  coefficients  in  Zones  I  and  II,  1/m 
x  -  cross -shore  coordinate  directed  positive  offshore,  m 
K  -  sand  transport  rate  coefficient,  m*/N 

D  -  wave  energy  dissipation  per  unit  water  volume,  N-m/m3-sec 
Deq  “  equilibrium  wave  energy  dissipation  per  unit  water  volume, 
N-m/m3-sec 

e  -  slope -related  sand  transport  rate  coefficient,  m2/sec 
h  -  still -water  depth,  m 

The  subscripts  b  ,  p  ,  z  ,  and  r  stand  for  quantities  evaluated  at  the 
break  point,  plunge  point,  end  of  the  surf  zone,  and  runup  limit,  respec¬ 
tively.  Different  spatial  decay  coefficients  are  used  in  Zones  I  and  II, 
denoted  by  the  subscripts  1  and  2,  to  describe  the  decrease  in  sand  transport 
rate  with  distance. 

52.  Empirical  expressions  for  the  spatial  decay  coefficients  were 
derived  from  large  wave  tank  data  (Report  1) .  The  decay  coefficient  in  Zone  I 
was  empirically  related  to  median  grain  size  (D50)  and  breaking  wave  height 
(Hjj)  as  (with  Dgg  in  mm  and  Hjj  in  m) 


(6) 


For  Zone  II,  limited  data  suggested  that 

\z  -  0.2  \  (7) 

In  calculating  the  transport  rate  in  Zones  I  and  II,  the  transport  rate  is 
first  determined  at  the  plunge  point  from  Equation  4,  and  then  the  exponential 
decay  rates  are  applied  seaward  in  the  respective  zones. 
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53.  The  cross -shore  transport  rate  in  zones  of  fully  broken  waves  is 
related  to  the  wave  energy  dissipation  per  unit  water  volume.  This  type  of 
transport  formula  has  previously  been  used  by  Moore  (1982),  Kriebel  (1982, 
1986),  and  Kriebel  and  Dean  (1985).  Analysis  of  the  large  wave  tank  data 
substantiated  this  type  of  transport  relationship  in  zones  of  broken  waves 
(Larson  1988,  Larson  and  Kraus  1989a).  The  transport  relationships  in  the 
other  zones  are  empirical  and  based  directly  on  the  data  from  the  wave  tank 
experiments.  For  example,  in  the  swash  zone  the  transport  rate  is  designed  to 
decrease  linearly  from  the  end  of  the  surf  zone  to  the  runup  limit.  This 
property  of  the  transport  rate  distribution  was  manifested  in  many  large  wave 
tank  cases,  occurring  with  both  onshore  and  offshore  transport.  Similar 
characteristics  have  been  noted  on  a  natural  beach  where  the  foreshore  changed 
uniformly  during  erosional  and  accretionary  events  (Seymour  1987).  However, 
no  direct  measurements  have  been  made  of  transport  rates  on  the  foreshore,  and 
transport  in  this  region  of  the  profile  deserves  further  study. 
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PART  III:  WAVE  MODEL 


54.  The  wave  model  provides  input  to  calculate  cross -shore  sand  trans¬ 
port  rates  from  which  profile  change  is  obtained.  In  this  chapter,  an 
overview  is  given  of  the  capabilities  and  limitations  of  the  wave  model 
including  a  summary  of  the  governing  equations.  The  equations  are  expressed 
in  an  explicit  finite-difference  form  and  presented  together  with  the  numer¬ 
ical  solution  method.  Finally,  the  interaction  of  the  wave  model  with  other 
modules  in  the  full  profile  change  model  is  discussed  with  reference  to  input 
and  output  requirements. 

55.  The  wave  model  is  a  generalized  version  of  the  model  described  in 
Report  1.  Significant  improvements  include: 

&.  Use  of  complete  linear  wave  theory  expressions  rather  than 
shallow-water  wave  approximations. 

b.  Addition  of  an  option  for  specifying  an  arbitrary  incident  wave 
angle  with  respect  to  the  bottom  contours.  The  bottom  contours 
are  assumed  to  be  straight  and  parallel. 

£.  Use  of  a  more  general  numerical  solution  scheme  based  on  the 
energy  flux  equation. 


Mofel  -S.apflbiUtje?  end  Limitations 


56.  Linear  wave  theory  is  used  to  determine  wave  characteristics 
across -shore  from  deep  water  or  a  specified  water  depth  offshore  to  the  break 
point.  Shoreward  of  the  break  point,  a  slightly  generalized  form  of  the  wave 
decay  model  of  Dally,  Dean,  and  Dalrymple  (1984,  1985)  calculates  the  wave 
height  distribution.  The  point  of  incipient  wave  breaking  is  determined  from 
an  empirical  criterion  expressed  in  terms  of  the  surf  similarity  parameter 
defined  with  the  deepwater  wave  steepness  and  the  local  beach  slope  seaward  of 
the  breakpoint,  evaluated  over  one -third  of  the  wavelength. 

57.  An  input  wave  height,  period,  and  incident  wave  angle  must  be  given 
either  in  deep  water  or  at  a  specified  depth.  The  model  uses  mean  wave  height 
to  calculate  transport  direction  and  significant  wave  height  to  calculate  the 
break  point  and  transport  rate  (Report  1).  Wave  information  is  transformed  by 
shoaling  and  refraction  over  straight  and  parallel  bottom  contours  to  the  most 
seaward  point  of  the  model  grid.  Energy  dissipation  in  the  viscous  layer 
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along  the  bottom  is  neglected  since  the  effect  on  wave  height  is  negligible  in 
comparison  to  the  amount  of  energy  dissipation  in  the  surf  zone. 

58.  Calculation  of  the  wave  height  across-shore  starts  in  the  most 
seaward  calculation  cell  and  proceeds  shoreward  cell  by  cell.  The  local  wave 
height  and  angle  with  respect  to  the  bottom  contours  are  determined  from 
conservation  of  energy  flux  and  Snell's  Law.  In  regions  of  wave  breaking, 
energy  dissipation  due  to  increased  turbulence  in  the  water  column  is  included 
in  the  energy  flux  equation.  It  is  possible  to  specify  bottom  contours  with 
an  orientation  that  is  not  parallel  to  the  baseline  (y-axis) . 

59.  Energy  dissipation  in  the  surf  zone  is  assumed  proportional  to  the 
excess  energy  flux  beyond  a  stable  energy  flux  that  is  a  function  of  water 
depth  (Horikawa  and  Kuo  1966).  Wave  properties  in  the  surf -zone  are  calcu¬ 
lated  using  linear  wave  theory  without  invoking  shallow-water  approximations. 
The  wave  height  calculation  is  carried  out  to  a  water  depth  that  is  less  than 
the  depth  where  the  surf  zone  ends  in  accordance  with  the  definition  of  the 
different  hydrodynamic  regions  in  the  nearshore.  The  depth  at  the  shoreward 
end  of  the  surf  zone  is  roughly  equal  to  the  seaward  limit  of  the  backrush. 

60.  The  wave  model  allows  wave  reformation  to  take  place  and  the  possi¬ 
bility  of  multiple  break  points  across  the  nearshore  zone.  Wave  reformation 
occurs  if  the  local  wave  height  reaches  the  stable  wave  height,  implying  no 
wave  energy  dissipation.  After  energy  dissipation  ceases,  shoaling  and 
refraction  become  dominant  until  wave  breaking  occurs  and  dissipation  is 
initiated  once  more. 

61.  Generalizing  the  breaker  decay  model  of  Dally,  Dean,  and  Dalrymple 
(1984,  1985),  the  two-dimensional  equation  for  conservation  of  energy  flux 
incorporating  energy  dissipation  associated  with  wave  breaking  may  be  written 


d  d 

—  (F  cose)  +  —  (F  sine)  -  ~  (F  -  Fs) 
8x  oy  d  s 


(8) 


where 

x  -  cross-shore  coordinate,  positive  directed  seaward,  m 
F  -  wave  energy  flux,  N-m/m-sec 
6  -  wave  angle  with  respect  to  the  bottom  contours 
y  «  longshore  coordinate,  m 
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k  -  empirical  wave  decay  coefficient 
d  -  total  water  depth,  m 
Fs  -  stable  wave  energy  flux,  N-m/m-sec 
The  total  water  depth  d  is  given  by 

d  -  h  +  r\  (9) 

where  t|  is  the  mean  water  surface  elevation  (setup  or  setdown)  produced  by 
wave  motion.  The  wave  energy  flux  is  given  by 

F  -  ECg  (10) 


where 

E  -  wave  energy  density,  N-m/m2 
Cg  -  wave  group  speed,  m/sec 

The  wave  energy  density  is  written  using  linear  wave  theory  as 

E  -  i  pgH2  (11) 

where 

p  -  density  of  water,  kg/m3 
H  -  wave  height,  m 

Assuming  the  wave  conditions  to  be  uniform  alongshore  and  the  bottom  contours 
to  be  straight  and  parallel,  Equation  8  may  be  reduced  to 

gj  (F  cos0)  -  g  (F  -  Fs)  (12) 

62.  Wave  group  speed  is  related  to  wave  phase  speed  C  through  a 
factor  n  which  is  a  function  of  the  water  depth  and  the  wavelength  L  (or 
wave  period  T  )  as 


C 


g 


-  nC 


(13) 
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where 


(14) 


The  wave  phase  speed  is  determined  through  the  dispersion  relationship 


C  -  CQ  tanh 

where 


(15) 


(16) 


in  which  CQ  is  the  wave  phase  speed  in  deep  water.  Equation  14  is  an 
implicit  equation  that  must  be  solved  numerically  or  by  rational  approxima¬ 
tions  since  C  contains  a  dependence  on  the  wavelength,  that  is,  C  -  L/T  . 

63.  The  wave  decay  coefficient  controls  the  rate  of  energy  dissipation, 
whereas  the  stable  energy  flux  determines  the  amount  of  energy  dissipation 
necessary  for  stable  conditions  to  occur  once  breaking  is  initiated.  In  this 
respect,  the  stable  wave  condition  refers  to  a  state  in  which  energy  dissipa¬ 
tion  during  breaking  ceases,  allowing  waves  to  reform.  The  stable  energy  flux 
may  be  expressed  as 


Fs  -  EsCg  (17) 

where  Es  -  stable  wave  energy  density,  N-m/m2  .  The  stable  wave  energy  flux 
corresponds  to  a  stable  wave  height  that  is  a  function  of  the  water  depth 

Hs  =  Td  (18) 

where  Hs  is  the  stable  wave  height  and  T  is  the  stable  wave  height 
coefficient  (empirical) . 
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64.  Two  empirical  coefficients  enter  in  the  breaker  decay  model,  k 
and  T  .  However,  the  values  of  these  coefficients  seem  to  vary  little 
over  a  wide  range  of  conditions,  both  in  the  laboratory  and  on  natural 
beaches.  On  the  basis  of  small-  and  large-scale  tank  data,  Dally,  Dean,  and 
Dalrymple  (1984)  recommended  the  values  of  k  -  0.15  and  T  -  0.40  for  use 
as  averages  for  a  varying  bottom  slope.  Ebersole  (1987)  evaluated  the 
performance  of  the  breaker  decay  model  using  field  measurements  and  found  good 
agreement  for  engineering  applications. 

65.  Report  1  verified  the  applicability  of  the  breaker  decay  model  for 
a  barred  beach  profile  using  the  large -wave  tank  measurements  of  Kajima  et  al. 
(1983) .  Values  obtained  for  the  two  empirical  coefficients  were  in  general 
agreement  with  values  obtained  by  the  aforementioned  authors.  An  apparent 
dependence  on  beach  slope  for  the  coefficients  was  noted,  but  no  specific 
relationships  could  be  developed. 

66.  Kraus  and  Larson  (in  preparation)  evaluated  the  performance  of  the 
breaker  decay  model  for  waves  approaching  the  beach  at  an  angle  with  respect 
to  the  bottom  contours  for  both  laboratory  (Visser  1982)  and  field  data  (Wu, 
Thornton,  and  Guza  1985).  In  these  cases,  the  coefficients  were  fixed  at  the 
standard  values  of  k  -  0.15  and  T  -  0.40.  In  all  cases,  the  breaker  decay 
model  produced  satisfactory  agreement.  For  the  field  data,  a  statistical 
approach  was  used  in  which  wave  heights  were  randomly  selected  from  a  Rayleigh 
distribution  for  a  large  number  of  waves,  and  the  calculated  wave  height 
distributions  were  averaged  to  obtain  a  representative  distribution. 

67.  As  waves  propagate  toward  shore,  a  flux  of  momentum  (radiation 
stress)  arises  and  causes  displacement  of  the  mean  water  elevation  if  changes 
in  this  flux  occur  due  to  shoaling  or  breaking.  Seaward  of  the  break  point, 
shoaling  produces  an  increase  in  wave  height,  which  causes  a  corresponding 
increase  in  the  momentum  flux.  This  flux  increase  is  balanced  by  lowering  of 
the  mean  water  elevation,  called  setdown  (Longuet-Higgins  and  Stewart  1963). 
Inside  the  surf  zone,  as  waves  continue  to  break  and  decrease  in  height,  the 
momentum  flux  decreases ,  and  an  increase  In  mean  water  elevation  occurs 
(setup) . 

68.  Displacement  of  the  mean  water  surface  (setup  or  setdown)  may  be 
determined  from  the  momentum  equation 
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(19) 


dsxx 

dx 


where  Sxx  is  the  radiation  stress  component  directed  onshore  (N/m) .  The 
radiation  stress  component  Sxx  is  given  by  linear  wave  theory  for  an 
arbitrary  wave  angle  of  incidence  as 


£n(ec 


Sxx  “  r  PS1*2  |n(eos29  +  1) 

O 


(20) 


69.  The  energy  flux  conservation  equation  (i  .aation  8),  the  wave 
dispersion  relationship  (Equation  15) ,  and  the  cross-shore  momentum  equation 
(Equation  19)  are  not  sufficient  to  specify  the  wave  height  across-shore  if 
the  waves  arrive  at  an  oblique  angle  to  the  bottom  contours.  In  this  case,  an 
additional  equation  must  be  solved  that  expresses  the  conservation  of  wave 
number.  The  wave  number  conservation  equation  is  written  in  the  general  form, 
not  taking  into  account  diffraction,  as 


3  fsinOl  3  fcosel 

k  rd  '  *y  nd  ~ 


(21) 


70.  For  straight  and  parallel  bottom  contours,  implying  no  alongshore 
variation  in  wave  properties,  Equation  21  may  be  simplified  to  yield 


(22) 


Integration  of  Equation  22  gives  Snell's  Law,  which  states  that  the  quantity 
sin0/L  is  constant  across-shore.  Equations  8,  15,  20,  and  22  allow  the  wave 
height  to  be  determined  once  boundary  conditions  have  been  established. 

71.  As  mentioned  previously,  the  wave  height,  period,  and  incident 
angle  are  calculated  at  the  seaward  end  of  the  profile  from  deepwater  or 
equivalent  conditions,  providing  the  necessary  boundary  conditions.  The  end 
of  the  grid  should  be  located  seaward  of  the  break  point  for  all  simulations. 
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Fiom  Che  wave  properties,  it  is  possible  to  determine  analytically  the  setdown 
at  the  most  seaward  point  (Longuet-Higgins  and  Stewart  1962). 


11  - 


nH2 


4L  sinh 


(23) 


72.  Also,  a  condition  for  the  initiation  of  breaking  is  needed, 
indicating  the  location  at  which  energy  dissipation  due  to  breaking  begins  and 
the  wave  decay  coefficient  becomes  different  from  zero.  Breaking  occurs  when 
the  ratio  between  the  wave  height  and  the  water  depth  exceeds  a  certain  value 
which  is  a  function  of  the  surf  similarity  parameter  expressed  in  terms  of  the 
deepwater  wave  steepness  and  profile  slope  prior  to  breaking. 

7  -  1.14  50*21  (24) 


where 

7  -  breaker  ratio  -  H^/h^ 

6  -  surf  similarity  parameter 
The  surf  similarity  parameter  is  given  by 


f H  7  *  0  •  5 

5  -  tanp  Ml  (25) 

where  tanfi  is  the  local  beach  slope  seaward  of  the  break  point. 

73.  Equation  24  was  derived  in  Report  1  from  data  given  by  Kajima  et 
al.  (1983),  but  very  similar  relationships  were  previously  reported  in 
Singamsetti  and  Wind  (1980)  and  Sunamura  (1980b) .  Report  1  used  data  from 
prototype- scale  experiments,  whereas  the  two  latter  studies  were  based  on 
small-scale  laboratory  experiments.  Thus,  in  Equation  24  scale  effects  are 
absent,  although  it  should  be  observed  that  only  monochromatic  waves  were  used 
in  the  tests. 

74.  The  laboratory  data  of  Kajima  et  al.  (1983)  included  cases  with  a 
maximum  wave  period  of  12  sec,  although  most  cases  had  periods  in  the  range  of 
3  to  9  sec.  The  limited  data  available  for  longer  period  waves  indicate  that 
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Equation  24  should  be  used  with  care.  Also,  since  the  relationship  was 
determined  mainly  from  barred  profiles  exposed  to  monochr  matic  waves,  the 
profile  slope  seaward  of  the  break  point  was  steeper  than  what  is  typically 
found  in  the  field.  Thus,  Equation  24  tends  to  underestimate  the  breaker 
ratio  for  gentle  slopes  in  field  applications. 

75.  An  expression  for  energy  dissipation  due  to  wave  breaking  is  needed 
in  the  prediction  of  the  cross -shore  sand  transport  rate.  In  the  process  of 
wave  height  decay  after  breaking,  turbulent  motion  is  induced  in  the  surf 
zone.  The  turbulent  kinetic  energy  is  converted  into  heat  through  transfer  of 
energy  from  larger  to  smaller  eddies.  In  the  macroscale  approach  of  the  beach 
profile  model,  it  is  satisfactory  to  derive  energy  dissipation  from  wave 
height  decay  without  resolving  the  internal  structure  of  this  process. 

However,  the  difference  between  the  production  of  turbulent  kinetic  energy  and 
the  corresponding  dissipation  with  reference  to  the  time  scale  (Roelvink  and 
Stive  1989)  is  observed  as  well  as  the  process  of  energy  reordering  (Svendsen, 
Madsen,  and  Buhr  Hansen  1978) . 

76.  The  wave  energy  dissipation  D  per  unit  water  volume  is 

1  d 

D  -  —  (F  cose)  (26) 

d  dx 

Using  Equation  12,  the  following  relationship  for  D  is  obtained 

D  -  ~z  (F  -  Fg)  (27) 


Once  the  cross-shore  wave  height  distribution  is  calculated,  the  energy 
dissipation  per  unit  water  volume  in  the  surf  zone  is  determined  explicitly 
from  Equation  27.  Outside  the  surf  zone,  k  is  set  to  zero,  and  no  energy 
dissipation  takes  place  since  bottom  friction  is  neglected. 
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Numerical  Solution  Scheme 


77.  An  explicit  finite-difference  scheme  is  used  to  solve  the  equations 
for  determining  cross-shore  wave  height.  In  the  following  section,  an  outline 
is  given  of  the  numerical  solution  scheme  and  the  difference  equations 
employed.  Figure  5  is  a  definition  of  the  grid  and  illustrates  where  the 
various  quantities  in  the  difference  equations  are  defined. 

78.  In  the  model,  a  quasi -stationary  approach  is  applied,  implying  that 
a  steady- state  solution  of  the  wave  height  distribution  is  determined  at  every 
time  step.  Propagation  of  individual  waves  is  not  described,  but  it  is 
assumed  that  the  model  time  step  is  significantly  larger  than  the  wave  period. 

79.  For  irregular  wave  conditions,  the  wave  height  used  at  each  time 
step  may  be  regarded  as  a  statistical  measure  that  represents  the  average 
properties  of  the  breaking  waves  during  the  time  step.  This  wave  height  is 
different  from  the  wave  height  needed  to  reproduce  the  average  wave  height 
distribution  during  a  time  step,  which  includes  all  broken  and  unbroken  waves 
at  specific  locations.  The  selected  wave  height  parameter  should  represent 
the  average  location  of  the  surf  zone,  giving  a  clear  distinction  between 
regions  of  broken  and  unbroken  waves.  From  experience  with  field  simulations, 
significant  wave  height  Is  satisfactory  in  this  respect  (Report  1) . 


qi-1  q  6  q£  +  1  q  0+2 


Quantities  evaluated  at  q-points:  q,  H,  0,  L,  F,  Cg 
Quantities  evaluated  at  h-points:  h,  tj,  Fg,  D 

Figure  5.  Definition  sketch  of  numerical  grid 
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80.  Numerical  calculations  start  at  the  seaward  end  of  the  grid  and 
proceed  onshore  through  an  explicit  solution  scheme  in  which  quantities  known 
at  a  specific  grid  point  are  used  to  determine  corresponding  quantities  at  the 
next  grid  point.  In  the  difference  equations  presented,  the  index  i  is  used 
to  denote  the  number  of  a  specific  grid  point.  A  staggered  grid  is  used  with 
different  quantities  taken  at  points  in  the  middle  of  calculation  cells  or  at 
the  boundaries  between  cells  (Figure  5).  The  primary  quantity  in  the  middle 
of  a  cell  is  water  depth,  and  these  grid  points  are  termed  h-points.  At  the 
boundaries  of  calculation  cells,  the  main  quantity  is  the  cros-- shore  trans¬ 
port  rate,  and  quantities  defined  here  are  termed  q-points. 

81.  As  previously  mentioned,  the  wave  height,  period,  and  incident 
angle  must  be  known  at  the  most  seaward  grid  point  (given  or  calculated  from 
some  reference  depth  offshore)  prior  to  starting  the  calculation.  From  these 
wave  properties,  wave  setdown  is  determined  from  Equation  23.  Profile  depth 
is  given  in  the  middle  of  each  calculation  cell,  whereas  setdown  and  setup  are 
calculated  at  q-points.  Depths  at  the  boundaries  of  the  cells  are  given  by 
linear  interpolation.  Also,  energy  flux  and  radiation  stress  in  the  direction 
of  the  waves  are  determined  from  knowledge  of  the  depth  and  wave  quantities  at 
the  first  q-point. 

82.  Moving  from  a  specific  q-point  to  the  next,  wave  refraction  is 
first  determined  if  the  incident  wave  approaches  with  an  angle  to  the  bottom 
contours.  From  Equation  22,  the  angle  between  the  wave  crests  and  the  bottom 
contours  at  the  next  shoreward  q-point  is  given  by 


9^  -  arcsin 


sin9 


't+l 


i+1 


(28) 


Note  that  grid  point  numbering  increases  in  the  seaward  direction  since  the  x- 
axis  points  offshore,  but  the  calculation  proceeds  from  offshore  in  the  shore¬ 
ward  direction.  Water  depth  at  grid  point  i  in  Equation  28  is  corrected 
with  setdown  (setup)  using  the  value  at  grid  point  i+1  to  avoid  use  of  an 
implicit  or  iterative  scheme.  Wavelengths  are  calculated  using  the  dispersion 
relationship  (Equation  15)  solved  by  a  Pade  approximation  (Hunt  1979) . 

83.  The  next  step  in  the  calculation  is  to  determine  the  energy  flux 
and  thus  the  wave  height.  Equation  10  is  written  in  difference  form  as 
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Fi  - 


cosQjl  +  0.5ACi 


[F1+1  (cos01+1  -  0.5ACi)  +  Ac  FSJ 


(29) 


where 


A„  - 


kAx 


Ci  hi  +  %+i 


(30) 


and  Ax  -  length  step  (m) .  The  stable  wave  energy  flux  is  determined  from 


V,  -  -  Pg  [r  (b,  ♦  wf  C?‘  * 

8  L  J  2 


(31) 


An  average  value  for  the  wave  group  speed  is  used  because  this  quantity  is 
taken  at  q-points,  and  the  stable  wave  energy  flux  is  evaluated  at  h-points. 
Seaward  of  the  break  point,  k  is  set  to  zero,  indicating  that  ACi  is  also 
zero,  since  no  energy  dissipation  due  to  breaking  occurs. 

84.  After  the  energy  flux  has  been  calculated  at  a  specific  point,  the 
corresponding  wave  height  is  determined  from  Equations  10  and  11  as 


Hi  “ 


(32) 


Using  the  wave  height,  radiation  stress  is  calculated  from  Equation  20  and 
setdown  (setup)  is  given  from  Equation  19,  expressed  in  difference  form  as 


(sxx)i41  ‘  (sxx\ 

’ll  "  Ui+i  +  - — - 1 
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(33) 


85.  At  every  calculation  step,  a  check  is  made  to  determine  if  wave 
breaking  has  occurred  according  to  Equation  24.  Once  breaking  is  initiated, 
the  wave  decay  coefficient  is  given  the  value  k  -  0.15  ,  and  energy  dissipa¬ 

tion  takes  place.  Also,  if  breaking  has  occurred,  a  check  is  made  to  deter¬ 
mine  whether  wave  reformation  is  possible,  that  is,  whether  the  wave  energy 
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flux  is  lower  than  the  stable  energy  flux.  The  stable  wave  height  coefficient 
is  set  to  r  -  0.40  ,  and  if  there  is  wave  reformation,  k  is  given  the 

value  zero.  An  arbitrary  number  of  surf  zones  with  intermediate  zones  of  wave 
reformation  may  occur. 

86.  The  breaker  ratio  is  determined  from  Equation  24  for  the  most 
seaward  break  point.  If  multiple  breaking  occurs,  a  constant  breaker  ratio 

7  -  1.0  is  applied  at  the  more  shoreward  break  points  to  avoid  extreme  values 
that  may  occur  by  use  of  Equation  24  at  steep-sloped,  small  bars  created  in 
the  inner  surf  zone.  The  shoreward  boundary  of  the  calculation  grid  is 
located  at  a  predetermined  water  depth  approximately  corresponding  to  the 
point  of  intersection  between  an  average  linear  sloping  foreshore  and  an  equi¬ 
librium  profile  for  a  given  sediment  diameter  or  the  location  of  a  seawall. 

87.  Wave  height  calculations  are  not  carried  out  using  the  exact 
topography  of  the  profile;  instead,  a  smoothed  profile  is  used.  The  reason 
for  this  is  that  the  response  of  waves  to  bottom  topography  takes  place  over  a 
certain  distance,  suggesting  that  shoaling  and  refraction  do  not  occur  instan¬ 
taneously.  For  example,  a  very  small  hump  (a  rock)  on  the  bottom  would  not 
cause  a  wave  to  shoal,  which  the  equations  would  predict  if  no  profile 
smoothing  were  applied.  The  smoothed  profile  is  used  only  to  compute  wave 
properties;  the  actual  calculated  profile  is  maintained  to  compute  transport 
rates  and  beach  profile  change. 

88.  The  number  of  calculation  cells  over  which  the  smoothing  is 
performed  (a  moving  average  scheme)  is  determined  based  on  the  breaking  wave 
height  as  3Hj>  .  The  averaging  distance  3Hjj  is  somewhat  arbitrary  and  is 
based  on  numerous  comparisons  of  model  predictions  made  in  Report  1.  The 
distance  over  which  a  wave  changes  markedly  in  shoaling  to  the  break  point  is 
not  known;  Smith  and  Kraus  (in  preparation)  have  recently  suggested  use  of  the 
distance  L/3,  where  L  is  the  average  local  wavelength,  a  qualitative 
observation  made  from  their  physical  model  tests  of  wave  breaking  over  bars. 

A  simple  predictive  formula  is  used  to  estimate  the  breaking  wave  height  at 
each  time  step  prior  to  determining  the  wave  height  distribution  (Larson  and 
Kraus  1989a)  as 
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(34) 
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PART  IV:  NUMERICAL  MODEL  OF  STORM- INDUCED  BEACH  CHANGE 


89.  The  model  SBEACH  consists  of  three  calculation  modules  that  are 
executed  consecutively  at  each  time  step  in  a  simulation.  The  modules 
calculate  wave  height  across-shore,  net  cross-shore  sand  transport  rate,  and 
profile  change,  respectively.  Wave  height  is  used  to  calculate  the  transport 
rate  from  which  profile  change  is  calculated.  The  bathymetry  is  then  revised, 
and  these  modules  are  stepped  through  successively  to  calculate  profile  change 
for  the  simulation  period. 

90.  Part  III  described  the  theoretical  background  and  numerical 
solution  scheme  of  the  wave  module.  Relationships  used  to  determine  the  net 
cross-shore  sand  transport  rate  were  given  in  Part  II.  In  this  chapter,  an 
overview  of  the  capabilities  and  limitations  of  the  profile  change  model  is 
given  with  a  discussion  of  the  input  data  requirements  and  types  of  output 
produced  by  the  model.  Finally,  the  numerical  solution  scheme  for  calculating 
profile  change,  boundary  conditions,  model  stability,  and  typical  time 
interval  and  length  step  is  described. 

Model  Capabilities  and  Limitations 

91.  A  fundamental  assumption  of  the  numerical  model  is  that  breaking 
waves  are  the  major  cause  of  profile  change,  and  thus  the  major  part  of  cross¬ 
shore  sand  transport  takes  place  in  the  surf  zone.  The  volume  of  sand 
transported  is  determined  by  the  wave  energy  dissipation  per  unit  water  volume 
in  the  surf  zone  (dependent  upon  wave  conditions  and  profile  shape)  and  by 
empirical  functional  relationships  obtained  from  the  large  wave  tank  data  in 
other  zones  along  the  profile.  The  direction  of  transport  is  give  by  an 
empirical  predictive  criterion  (Equation  1) ,  and  calculation  of  the  fluid 
velocity  field  in  the  surf  zone  is  not  required. 

92.  Wave  height,  wave  period,  and  water  level  can  be  specified  at  each 
time  step  to  account  for  variable  forcing  conditions.  Water  level  change 
appears  in  the  computations  as  a  translation  of  the  origin  of  the  vertical 
coordinate  axis;  the  water  level  is  shifted  back  to  the  original  level  prior 
to  entering  the  mass  conservation  equation. 
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Q3.  The  numerical  model  incorporates  the  principle  of  superposition  in 
the  sense  that  random  waves,  as  exist  in  the  field,  are  represented  by 
consecutive  individual  waves  appearing  in  a  time  series  of  monochromatic 
waves.  This  assumption  is  appropriate  with  respect  to  the  engineering  utility 
of  the  model  where  a  macroscopic  quasi- stationary  approach  is  taken  in 
calculating  the  cross-shore  wave  height  distribution.  Also,  the  time  interval 
should  be  sufficiently  short  to  represent  the  many  different  waves  in  a  time 
series  of  random  waves. 

94.  For  the  large  wave  tank  data,  choice  of  input  wave  height  and 
period  for  the  computations  is  unambiguous  because  waves  of  constant  height 
and  period  were  used.  On  a  natural  beach,  the  forcing  parameters  show  apprec¬ 
iable  variation  in  time.  By  trial  calculations  using  field  data,  it  was 
concluded  that  significant  wave  height  gives  the  best  result  regarding  the 
transport  rate  and  corresponding  profile  change  as  compared  with  mean  or  root- 
mean-square  wave  height  (Report  1) .  An  intuitive  explanation  for  using 
significant  wave  height  is  that  energy  dissipation,  which  determines  the 
transport  rate,  increases  with  the  square  of  the  wave  height  implying  that 
higher  waves  give  a  relatively  larger  contribution.  Thus,  as  a  representative 
wave  height,  a  larger  measure  should  be  used  rather  than  the  mean  wave  height. 

95.  Mimura,  Otsuka,  Watanabe  (1986)  conducted  small-scale  laboratory 
experiments  using  irregular  waves  in  order  to  determine  the  representative 
monochromatic  wave  height  for  describing  various  cross -shore  transport 
phenomena.  They  found  that  the  overall  profile  response  as  erosional  or 
accretionary  was  best  described  by  the  mean  wave  height.  This  finding  is  in 
agreement  with  the  results  in  Report  1  as  displayed  by  the  criterion  separ¬ 
ating  erosion  and  accretion  (Equation  1  and  Figures  1  and  2) .  The  study  by 
Mimura,  Otsuka,  and  Watanabe  did  not  provide  clear  guidance  concerning  choice 
of  a  representative  wave  for  calculating  the  net  cross-shore  transport  rate, 
but  both  the  mean  and  significant  wave  height  appeared  to  be  applicable. 

Their  transport  relationship  (Watanabe  1982)  is  based  on  the  Shields  para¬ 
meter,  which  is  conceptually  and  functionally  different  from  the  wave  energy 
dissipation  approach  used  in  Report  1. 

96.  Median  grain  diameter  for  the  studied  profile  is  a  key  parameter  in 
calculating  the  net  cross-shore  transport  rate  and  resulting  profile  change. 
This  representative  sand  size  determines  the  equilibrium  energy  dissipation 
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which  appears  in  the  transport  rate  formula  in  the  zones  of  broken  waves 
(Equation  4) .  Moore  (1982)  presented  a  relationship  between  De q  and  D5Q 
in  graphical  form  based  on  field  profiles,  and  Dean  (1987,  1990)  re-expressed 
the  relationship  in  terms  of  fall  speed.  As  described  in  Report  1,  the 
numerical  model  gave  the  best  overall  fit  for  large  wave  tank  data  if  0.75  De q 
was  used  in  the  simulations  rather  than  Deq  .  This  proved  satisfactory  for 
field  data  as  well. 

97.  In  the  numerical  model  it  is  possible,  in  principle,  to  specify  a 
different  grain  size  and  thus  equilibrium  energy  dissipation  in  every  calcula¬ 
tion  cell.  However,  since  the  movement  of  individual  grains  is  not  tracked, 
the  median  grain  diameter  will  remain  the  same  throughout  the  simulation  in  a 
specific  cell,  which  is  unrealistic  considering  grain  movement.  Some  beach 
profiles  exhibit  grain  size  characteristics  with  distinctly  different  regions 
across -shore.  If,  for  instance,  the  dune  region  has  a  median  grain  diameter 
different  from  the  subaqueous  part  of  the  profile,  it  may  be  permissible  to 
schematically  represent  these  conditions.  In  this  case,  every  cell  located  in 
the  dune  region  is  assigned  one  grain  size,  whereas  the  cells  located  seaward 
are  given  a  different  grain  size.  Using  this  technique,  it  is  implicitly 
assumed  that  as  grains  pass  from  the  dune  region  to  the  seaward  region,  the 
transported  amount  exchanged  between  regions  is  so  small  that  it  does  not 
alter  the  original  distribution  of  grain  size. 

98.  Development  of  the  numerical  model  was  directed  towards  predicting 
profile  change  during  erosional  and  accretionary  events.  Two  main  areas  of 
application  of  the  model  are  estimation  of  profile  erosion  and  shoreline 
retreat  during  storm  or  hurricane  events  and  initial  adjustment  of  beach  fills 
to  changes  in  water  level  and  waves.  Useful  descriptive  quantities  are  eroded 
volume  above  a  certain  contour  and  retreat  distance  of  that  contour.  The 
program  calculates  these  quantities  for  an  arbitrarily  specified  contour  and 
selected  time  step.  These  applications,  together  with  the  properties  of  the 
large  wave  tank  experiments,  emphasized  model  development  focusing  on 
erosional  processes. 
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Profile  Change  Model 


99.  Profile  change  is  calculated  from  the  mass  conservation  equation 
using  the  net  transport  rate  distribution  determined  from  Equations  2  through 
5.  If  sand  sources  or  sinks  are  absent,  the  mas?  conservation  equation  is 


3h  3q 
3t  "  3x 


(35) 


The  boundary  conditions  are  zero  sand  transport  rate  at  the  runup  limit  and  at 
the  seaward  depth  signifying  the  onset  of  significant  sand  movement.  The 
seaward  boundary  of  the  calculation  grid  (limiting  depth  for  sand  motion)  is 
effectively  determined  by  the  exponential  decay  in  Zone  I,  taken  to  be  a  point 
where  the  transport  rate  has  decreased  to  a  small  value. 

100.  As  the  wave  height  distribution  is  calculated  across -shore  for  a 
given  time  step,  the  location  of  the  boundaries  between  the  different  sand 
transport  zones  is  determined.  The  location  of  the  break  point  is  given 
directly  from  the  wave  height  calculations,  and  the  plunge  point  is  computed 
as  3Hb  based  on  findings  by  Galvin  (1969)  and  Svendsen  (1987).  The  landward 
boundary  of  the  surf  zone  is  defined  as  the  point  of  tangency  between  a 
linearly  sloping  beach  face  and  an  equilibrium  profile  for  a  given  median 
grain  diameter,  whereas  the  runup  limit  is  given  by  a  formula  for  predicting 
the  active  profile  height  developed  from  the  large  wave  tank  data  (Report  1) . 
The  expression  for  active  profile  height  is  similar  to  the  formula  provided  by 
Hunt  (1959),  where  the  runup  height  normalized  with  the  incident  wave  height 
is  a  function  of  the  surf  similarity  parameter. 

101.  The  active  subaerial  profile  height  ZR  is  given  by 

ZR  -  1.47  50-79  (36) 

The  surf  similarity  parameter  F  is  defined  in  Equation  25,  and  the  slope  to 
be  used  in  Equation  36  is  the  average  beach  slope  in  the  surf  zone.  Origi¬ 
nally  Equation  36  was  based  on  the  initial  slope  in  the  large  wave  tank  tests. , 
but  application  of  the  model  to  field  data  showed  that  using  the  average  surf 
zone  slope  gave  satisfactory  results.  The  runup  height  has  little  influence 


42 


on  bar  size  and  movement,  but  it  does  define  the  shoreward  limit  for  cross¬ 
shore  sand  movement  (Report  1) .  The  amount  of  erosion  on  the  foreshore  is 
similar  even  for  significant  variations  in  the  runup  height,  but  the  shape  of 
the  foreshore  is  quite  different.  Transport  at  the  landward  end  of  the  surf 
zone  determines  how  much  material  is  transported  from  or  to  the  foreshore. 
Higher  runup  implies  a  larger  area  from  which  material  can  be  taken  and  thus  a 
different  resulting  foreshore  shape  compared  with  that  given  by  a  smaller 
runup  height. 

102.  From  knowledge  of  the  locations  defining  the  various  transport 
zones,  the  transport  rate  distribution  can  be  calculated  from  Equations  2 
through  5.  Several  break  points  may  occur  along  a  profile  if  wave  reformation 
takes  place,  leading  to  several  zones  of  Type  II  and  III  transport.  To 
determine  the  transport  rate  distribution,  sand  transport  is  first  calculated 
in  zones  of  fully  broken  waves  (Zone  III)  according  to  Equation  4,  written  in 
difference  form 


qi  -  k 


j~(Qi  +  Pi-i) 


"  De<?  +  KAx  (Hi  ’  hi_l) 


(37) 


The  definition  of  the  calculation  grid  is  given  in  Figure  5.  Energy  dissipa¬ 
tion  per  unit  water  volume  is  taken  at  h-points  (in  the  middle  of  a  calcula¬ 
tion  cell)  making  it  necessary  to  average  D  when  calculating  the  transport 
rate. 

103.  Calculation  of  the  transport  rate  in  zones  of  fully  broken  waves 
defines  the  boundaries  of  Zone  III,  from  which  the  transport  rate  may  be 
calculated  in  the  other  zones.  After  these  values  are  determined  at  the 
plunge  point  and  at  the  end  of  the  surf  zone,  Equations  2,  3,  and  5  are 
applied  to  completely  define  the  transport  rate  distribution.  In  discretizing 
the  mass  conservation  equation,  transport  rate  distributions  from  two  time 
levels  are  used.  The  superscript  k  is  used  to  denote  a  specific  time  level, 
and  the  difference  equation  is  expressed  as 


At 


1  fffi  -  <6n 

2  [  Ax 


(38) 
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where  At  is  the  time  step  (sec).  Since  the  transport  rate  distribution  is 
determined  from  different  relationships  depending  on  the  zone  of  transport, 
its  spatial  derivative  will  generally  be  discontinuous  at  the  boundaries 
between  the  zones.  To  obtain  a  smoother,  more  realistic  transport  distribu¬ 
tion,  a  3 -point  filter  is  applied  to  the  calculated  transport  rates  on  grid 
cells  away  from  boundaries  by 

-  0.25  qL_i  +  0.50  q£  +  0.25  qi+1  (39) 

where  the  prime  (')  denotes  the  smoothed  transport  rate. 

104.  In  some  large  wave  tank  cases,  the  local  beach  slope  exceeded  the 
angle  of  repose  due  to  constant  waves  and  water  level.  Monochromatic  waves 
caused  erosion  and  accumulation  to  occur  at  the  same  locations  along  the 
profile  for  a  sufficiently  long  duration  to  produce  excessive  profile  steepen¬ 
ing.  Steepest  profile  slopes  occurred  during  erosional  conditions  on  the 
shoreward  side  of  the  bar  just  seaward  of  the  berm  crest.  Local  profile 
slopes  exceeding  the  angle  of  repose  rarely  develop  on  a  natural  beach  because 
wave  height,  wave  period,  and  water  level  are  continuously  varying. 

105.  The  concept  of  avalanching,  as  discussed  by  Allen  (1970),  was 
incorporated  in  the  numerical  model,  and  a  routine  was  included  to  account  for 
transport  induced  by  slope  failure.  Allen  recognized  two  different  limiting 
profile  slopes,  the  angle  of  initial  yield  (angle  of  repose)  and  the  residual 
angle  after  shearing.  If  the  local  slope  exceeds  the  angle  of  initial  yield 
iji^y  ,  material  is  redistributed  along  the  slope,  and  a  new  stable  slope  is 
attained,  known  as  the  residual  angle  after  shearing  i^ra  •  From  analysis  of 
maximum  local  profile  slopes  from  the  large  wave  tank  data,  considerably  lower 
values  were  obtained  for  i^iy  and  t|»ra  than  reported  by  Allen.  From  his 
experiments  with  natural  sands  (diameters  ranging  from  0.27  to  3.17  mm),  he 

JL 

noted  an  angle  of  about  48  deg  caused  avalanching,  whereas  an  angle  of  33  deg 
seemed  to  be  the  stable  slope  after  avalanching  had  ceased. 

106.  Analysis  of  the  large  wave  tank  data  (Report  1)  showed  smaller 
values  of  the  maximum  slope,  inferred  to  be  the  angle  of  initial  yield,  with  a 


*  A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  5. 
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mean  value  of  28  deg.  The  average  value  for  the  residual  angle  after  shearing 
was  determined  to  be  22  deg;  however,  this  conclusion  is  not  firm  since  some 
profile  steepening  undoubtedly  occurred  between  the  time  of  avalanching  and 
the  next  profile  survey.  Allen  determined  from  experiments  that  the  differ¬ 
ence  between  i|>£y  and  i|>ra  (dilatation  angle)  is  in  the  range  of  10  to  15 
deg  for  sand.  Therefore,  in  the  numerical  model,  if  ij’iy  is  set  to  28  deg 
and  i|>ra  is  set  to  18  deg,  a  dilatation  angle  of  10  deg  is  used. 

107.  If  the  local  slope  exceeds  4*iy  in  the  model,  material  will  be 

redistributed  in  the  neighboring  cells  so  that  the  slope  everywhere  is  less 
than  or  equal  to  i|>ra  .  Computationally  this  is  performed  by  checking  the 
local  slope  iC  at  each  time  step  along  the  grid,  and,  if  i|>  >  *  avalanch¬ 

ing  is  initiated.  Two  conditions  determine  the  redistribution  of  material; 
mass  should  be  conserved  and  the  residual  angle  after  shearing  is  i|>ra  . 

Figure  6  is  a  definition  sketch  for  avalanching  where  the  local  slope  between 
cells  1  and  2  exceeds  ij’iy  •  Depths  before  avalanching  are  denoted  with  h  , 
whereas  new  depths  after  sand  redistribution  are  indicated  with  a  prime. 

108.  The  number  of  cells  N  over  which  avalanching  occurs  is  not  known 
a  priori  and  has  to  be  found  iteratively.  Initially,  it  is  assumed  that  the 
avalanching  takes  place  between  cells  1  and  2.  Material  is  moved  from  cell  1 
to  2  in  order  to  obtain  the  slope  t|>ra  between  the  cells,  and  the  new  slope 
between  cells  2  and  3  is  checked  to  see  if  it  exceeds  ipy^  .  If  this  is  the 
case,  the  calculation  is  repeated,  now  incorporating  cells  1  to  3,  using  the 
mass  conservation  equation  and  assigning  the  slope  i|>ra  to  all  slopes  between 
cells.  This  procedure  is  continued  until  the  slope  between  cells  N  and  N+l 
is  less  than  4>ra  . 

109.  For  the  general  case  where  sand  is  redistributed  in  N  calcula¬ 
tion  cells  included  in  avalanching,  the  mass  conservation  equation  is  written 

N 

l  Ah,  -  0  (40) 

where 

Ahi  -  hi  -  ^  (41) 
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Figure  6.  Definition  sketch  for  avalanching  calculation 


Here  the  cell  numbering,  denoted  by  i  ,  starts  from  the  cell  where  avalanch¬ 
ing  is  initiated  and  is  taken  in  the  direction  of  sand  redistribution.  Also, 
depth  changes  have  to  be  taken  with  the  proper  sign.  The  other  criterion 
necessary  to  obtain  the  solution  is  that  after  avalanching,  the  slope  between 
calculation  cells  must  be  ij>ra  .  If  the  length  of  the  calculation  cell  is 
Ax,  the  residual  angle  after  shearing  is  expressed  as 


^ra 


(42) 


Since  s|>ra  is  constant,  the  difference  in  depth  between  neighboring  cells  Ah 
has  to  be  a  constant  given  by 

Ah  -  *ra  Ax  (43) 

110.  The  requirement  on  the  difference  in  depth  between  cells,  as  given 
by  Equation  43,  may  be  written  for  N-l  pairs  of  cells  included  in  the 
avalanching  as 
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Ah  -  (hi  +  Ahi)  -  (h2  +  Ah2) 

Ah  -  (h2  +  Ah2)  -  <h3  +  Ah3) 

Ah  -  +  Ahi)  -  (hi+1  +  Ahi+1) 

Ah  -  (h„_i  +  Ah,,^)  -  (hM  +  AhH) 


(44) 


This  system  of  equations  may  be  rewritten  by  consecutively  eliminating  the 
first  part  of  the  right-hand  side  of  each  equation  by  adding  the  previous 
equation.  The  result  is 


Ah  -  (h2  +  Ahx) 
2Ah  -  (hi  +  Aht) 
iAh  -  (hi  +  Ahj) 
(N-l)Ah  -  (hi  +  A^) 


(h2  +  Ah2) 

(h3  +  Ah3) 
(hi+i  +  4hi+1) 
(hH  +  AhH) 


(45) 


111.  From  Equation  45,  the  unknown  depth  changes  Ahi  may  be  expli¬ 
citly  identified  as 


Ah2  -  (hi  +  Ah^ 
Ah3  -  (hx  +  Ahx) 
Ahi+1  -  (hj  +  Ahi) 
Ah„  -  (hi  +  Ahi) 


h2  -  Ah 
h3  -  2Ah 
hi+i  -  iAh 
hN  -  (N-l)Ah 


(46) 


From  Equation  46,  it  is  possible  to  compute  all  the  depth  changes  once  the 
depth  change  in  the  first  cell  is  known.  To  determine  the  depth  change  in  the 
first  cell,  the  mass  conservation  equation  (Equation  40)  is  used.  The  sum  of 
all  changes  in  depth  must  be  zero 


N 

l  Ahi 


i«l 


N 

Ah:  +  £  Ahi  "  0 

i-2 


(47) 
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Replacing  the  summation  over  the  depth  changes  from  cell  2  to  N  using 
Equation  46  yields 


N  KM 

Aht  +  l  <hx  +  Ahj)  -  l  hx  -  l  Ah  - 
1-2  1-2  1-2 


(48) 


112.  From  Equation  48,  it  is  possible  to  explicitly  solve  for  the  depth 
change  in  cell  1. 


(N-l)  1  "  1 

Ah,  -  -  - h,  +  -  T  h,  +  -  (N-l)  Ah 

1  N  N  i-2  2 


(49) 


Depth  changes  in  cells  2  to  N  are  then  given  by 

AhA  -  hx  +  Ahj  -  ^  -  (i-1)  Ah  (50) 

113.  The  time  interval  in  the  numerical  model  typically  varies  between 
5  and  20  min,  whereas  the  length  step  is  typically  1  to  5  m.  The  time 
interval  and  length  step  are  not  independent  because  they  both  govern  the 
stability  of  the  numerical  solution  scheme.  Generally,  a  shorter  length  step 
requires  a  smaller  time  interval.  An  explicit  criterion  for  the  relation 
between  time  and  length  step  for  numerical  stability  could  not  be  determined 
due  to  the  nonlinearity  of  the  equations  and  the  smoothing  of  the  transport 
rate.  However,  it  was  empirically  found  that  use  of  the  slope -dependent  term 
in  the  transport  equation  in  the  broken  wave  zone  (Equation  4)  greatly 
improves  stability  of  the  solution.  If  this  term  is  omitted,  numerical 
oscillations  may  originate  from  the  shoreward  side  of  the  bar  as  the  profile 
slope  steepens  under  constant  waves  and  water  level. 
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PART  V:  MODEL  TESTS 


114.  Although  the  numerical  model  successfully  reproduced  beach  profile 
change  generated  in  large  wave  tanks,  including  bar  formation,  bar  movement, 
and  foreshore  erosion,  the  ultimate  test  is  to  verify  model  performance  for 
field  beaches.  This  chapter  reviews  previous  field  simulations  (Report  1) 
that  concentrated  on  bar  formation  and  movement  and  discusses  recent  field 
simulations  where  the  objective  was  to  test  the  model  for  reproducing  fore¬ 
shore  erosion  as  well  as  accretion.  In  the  field  tests,  documented  storm 
erosion  events  from  both  the  US  east  and  west  coasts  were  simulated.  The 
chapter  begins  with  a  summary  of  sensitivity  analyses  performed  to  investigate 
predictions  of  bar  features  to  changes  in  model  parameters. 

Sensitivity  Analyses  and  Previous  Model  Tests 

115.  The  objective  of  a  sensitivity  analysis  is  to  establish  the 
influence  of  model  parameters  on  simulation  results.  Report  1  describes 
simulations  performed  to  quantify  the  effect  of  variations  in  model  parameters 
on  the  prediction.  The  investigation  focused  on  bar  properties  such  as 
volume,  maximum  height,  and  location  of  the  center  of  mass.  Parameters 
investigated  were  the  transport  rate  coefficients  K  and  e  ,  breaker  decay 
model  coefficients  k  and  T  ,  and  equilibrium  energy  dissipation  Deq  .  The 
influence  of  wave  height,  wave  period,  and  runup  height  were  also  explored. 

116.  Grain  size  is  a  critical  input  parameter  that  determines  the 
equilibrium  energy  dissipation  and  thus  the  amount  of  material  that  is  moved 
by  the  waves  in  the  course  of  reaching  an  equilibrium  beach  shape.  Figure  7 
displays  the  time  evolution  of  bar  volume  for  various  grain  sizes  in  a 
hypothetical  erosional  case.  The  equilibrium  beach  profile,  defined  as  the 
profile  shape  for  which  no  net  transport  occurs  and  bar  growth  ceases  for  a 
given  wave  condition,  is  approached  more  rapidly  for  coarser  grain  sizes. 

This  is  expected  since  less  material  has  to  be  redistributed  by  waves  for  a 
larger  value  of  Deq  .  Grain  size  is  one  of  the  few  model  parameters  that 
significantly  affects  the  movement  of  the  bar  center  of  mass.  Figure  7  is 
typical  of  the  results  from  the  sensitivity  analysis,  and  similar  curves  are 
presented  in  Report  1  for  all  model  parameters  discussed  here. 
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Figure  7.  Evolution  of  bar  volume  for  different  grain 
sizes  for  hypothetical  simulations 


117.  The  transport  rate  coefficient  K  mainly  governs  the  time 
response  of  the  profile  and  has  much  less  influence  on  bar  Vv'ume  and  eroded 
foreshore  sand  volume .  A  large  value  of  K  implies  rapid  profile  response 
with  rapid  attainment  of  equilibrium.  The  coefficient  e  ,  controlling  the 
strength  of  the  slope -dependent  term  in  Equation  4,  does  not  greatly  influence 
the  time  response,  particularly  at  the  initial  phase  of  a  simulation,  and 
mainly  contributes  to  reduce  the  depth  of  the  trough  landward  of  the  bar.  Bar 
volume,  however,  is  affected  by  a  change  in  e  ,  where  a  larger  value  implies 
a  larger  bar  volume. 

118.  The  empirical  coefficients  in  the  breaker  decay  model  have  a 
marked  influence  on  profile  evolution  with  regard  to  bar  volume.  A  decrease 
in  k  causes  an  increase  in  bar  volume  that  may  be  physically  explained  by 
the  resultant  larger  wave  energy  dissipation  in  the  inner  part  of  the  surf 
zone.  A  smaller  k -value  produces  more  gradual  energy  dissipation,  allowing 
waves  to  retain  more  energy  to  be  dissipated  farther  inshore.  A  large  value 
of  the  stable  wave  height  coefficient  T  implies  a  flatter  equilibrium  beach 
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profile,  requiring  more  material  to  be  moved  before  equilibrium  is  attained 
and  thus  more  foreshore  erosion. 

119.  The  effect  of  changing  wave  height  and  wave  period  is  obvious, 
where  an  increase  in  either  quantity  implies  greater  incident  wave  energy  and 
more  erosion.  Runup  height  has  little  influence  on  bar  development  since 
about  the  same  amount  of  sand  is  eroded  along  the  profile  independent  of  runup 
height.  However,  the  shape  of  the  foreshore  can  be  substantially  different 
depending  on  runup  height.  Smaller  runup  height  causes  erosion  to  occur  over 
a  narrower  portion  of  the  foreshore  in  comparison  to  larger  runup;  a  more 
pronounced  scarp  develops  with  a  smaller  runup  height. 

120.  In  Report  1,  the  numerical  model  was  used  to  simulate  profile 
change  measured  at  the  FRF,  located  on  a  barrier  island  beach  facing  the 
Atlantic  Ocean  at  Duck,  NC.  High-quality  beach  profile  surveys  from  the  dune 
to  the  depth  of  significant  sand  transport  have  been  made  at  approximately 
2-week  intervals  for  more  than  4  years.  Four  different  shore-normal  lines 
were  monitored  together  with  more  frequent  wave  and  water  level  measurements. 
Gage  measurements  of  wave  height  and  period  were  made  at  6 -hr  intervals,  and 
water  level  is  recorded  every  6  min  but  available  at  hourly  intervals.  A 
portion  of  the  profile  data  set  is  given  in  Howd  and  Birkemeier  (1987). 

121.  Five  beach  profile  sets  were  chosen  for  supplying  simulation  data 
for  model  verification.  Each  set  contained  two  profile  surveys  with  a  time 
history  of  waves  and  water  levels  between  surveys.  Erosion  predominated  for 
all  five  events,  including  significant  offshore  movement  of  one  or  two  bars, 
but  shoreline  position  showed  little  or  no  movement.  Shoreline  stability  is 
characteristic  at  the  FRF  beach  and  may  be  caused  in  part  by  an  abundance  of 
coarse  material  on  the  foreshore.  The  data  set  from  Duck  was  thus  unsuitable 
for  testing  predictability  of  the  numerical  model  regarding  shoreline  retreat 
and  subaerial  beach  erosion,  but  the  data  were  considered  appropriate  for 
evaluating  bar  movement. 

122.  Potential  three-dimensionality  of  changes  on  a  natural  beach  must 
be  considered  in  simulating  profile  change.  Even  though  simulation  time 
periods  were  chosen  such  that  neighboring  profile  surveys  most  nearly  dis¬ 
played  similar  development,  some  three-dimensional  effects  were  present.  The 
degree  of  three-dimensionality  was  quantified  by  applying  the  mass  conserva¬ 
tion  equation  to  calculate  the  loss  or  gain  of  sand  between  consecutive 


51 


profiles.  A  survey  line  located  southeast  of  the  FRF  pier  was  used  in  the 
simulation  because  it  is  beyond  major  influence  of  the  pier,  and  the  bottom 
contours  are  essentially  straight  and  parallel  (Howd  and  Birkemeier  1987) . 

123.  Four  of  the  five  events  were  used  for  calibration,  where  the 
calibration  criterion  comprised  minimizing  the  sum  squares  of  the  difference 
between  measured  and  calculated  profile  change.  The  fifth  event  was  used  for 
model  verification  (Report  1) .  The  energy-based  significant  wave  height  (H,^ 
from  gage  620  located  in  18  m  of  water)  was  used  to  calculate  the  transport 
rate  distribution  together  with  the  peak  period  of  the  spectra.  Since  the 
measurement  intervals  of  waves  and  water  level  were  longer  than  the  model  time 
step,  required  intermediate  values  were  obtained  by  cubic  spline  interpola¬ 
tion.  A  typical  simulation  time  step  was  20  min,  and  the  length  step  was  5  m. 

124.  Essentially  the  same  parameter  values  obtained  from  the  large  wave 
tank  simulations  were  used,  although  some  modifications  had  to  be  made  for 
field  simulations.  The  most  important  difference  was  that  the  transport 
coefficient  K  had  to  be  decreased  to  obtain  optimal  agreement  for  field 
simulations.  The  optimal  value  of  K  for  field  simulations  was  about  half 
the  value  for  the  large  wave  tank  data.  Two  different  median  grain  diameters 
were  used  to  represent  sediment  characteristics  of  the  FRF  beach  along 
transect  188  (Howd  and  Birkemeier  1987).  A  median  grain  size  of  2.0  mm  was 
used  on  the  foreshore,  and  0.15  mm  was  used  for  points  seaward  of  the 
foreshore . 

125.  Figure  8a  shows  the  measured  wave  height,  wave  period,  and  water 
level  for  the  simulation  period,  developed  as  a  spline  fit  through  measured 
values  at  1-hr  intervals.  Figure  8b  illustrates  one  of  the  calibration  runs 
for  3-6  April  1984  when  a  single  bar  was  present  that  moved  about  60  m 
offshore  during  a  minor  storm  event.  Measured  initial  and  final  profiles  are 
given  together  with  the  final  calculated  profile.  Bar  location  was  rather 
well  predicted,  whereas  bar  volume  was  underestimated  and  the  trough  was  not 
sufficiently  pronounced.  Also  some  shoreline  retreat  was  calculated  in  spite 
of  the  larger  grain  size  on  the  foreshore. 

126.  Model  verification  gave  reasonable  agreement  between  measured  and 
calculated  profiles,  although  the  selected  event  encompassed  profiles  exhibit¬ 
ing  two  bar  features,  which  was  not  the  case  in  the  calibration  events.  As 
with  the  calibration  cases,  direction  of  bar  movement  and  location  were 


52 


Figure  8.  Calibration  of  numerical  model  against  field  data 

from  Duck,  NC 
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correctly  predicted,  whereas  the  size  of  the  bars  was  smaller  than  measured 
(Figure  9) .  S5.mulations  performed  subsequent  to  these  early  field  tests 
indicate  that  a  probable  reason  for  the  smaller  calculated  bar  size  is 
insufficient  variability  in  the  imposed  wave  input.  Greater  variability  would 
tend  to  shift  the  break  point  back  and  forth  to  produce  a  wider  bar.  Some 
variability  in  wave  and  water  level  conditions  is  reproduced  through  use  of  a 
time  series  of  different  regular  waves,  but  it  is  believed  that  introduction 
of  random  fluctuations  around  these  values  would  give  more  realistic  predic¬ 
tions.  The  absence  of  randomness  in  the  wave  description  is  probably  one  of 
the  factors  requiring  use  of  a  lower  transport  rate  coefficient  in  field 
simulations  as  compared  with  large  wave  tank  simulations.  A  random  variation 
in  wave  height,  for  example,  would  reduce  the  tendency  for  a  steady  trend  in 
profile  evolution,  thereby  requiring  a  larger  value  of  the  transport  rate 
coefficient  K  to  achieve  optimal  agreement. 

127.  Summarizing  experience  from  modeling  beach  profile  evolution  at 
Duck,  prediction  of  the  direction  of  bar  movement  and  bar  location  was 
successful,  although  calculated  bar  size  showed  less  agreement  with  observed 
profile  development.  Profile  response  was  best  reproduced  using  as  much 
information  as  possible  on  wave  height,  wave  period,  and  water  level  (as 
opposed  to  using  average  values  for  the  simulation  period),  lending  credibi¬ 
lity  to  the  predictive  capabilities  of  the  numerical  model.  The  success  of 
the  model  in  predicting  observed  bar  movement  supports  the  hypothesis  that 
breaking  waves  are  the  main  agent  for  cross -shore  sand  movement,  at  least  in 
the  outer  surf  zone.  The  steep  bar  and  foreshore  slopes  appearing  in  the 
large  wave  tank  simulations  were  not  encountered  in  field  simulations  because 
measured  changing  forcing  conditions  were  used.  Reproduction  of  milder  bar 
features  validates  the  concept  of  simulating  profile  change  by  random  waves 
through  a  summation  of  transport  events  produced  by  different  regular  waves 
acting  over  short  time  intervals. 
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Depth  (m)  Wave  Height  and  Water  Level  (m) 


a.  Variation  of  wave  height,  wave  period,  and  water  level 


Figure  9.  Verification  of  numerical  model  against  field  data 

from  Duck,  NC 
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Storm  Event  on  the  US  East  Coast 


128.  Several  profile  lines,  located  north  and  south  of  Manasquan  Inlet, 
NJ,  were  surveyed  by  the  US  Army  Engineer  District  (USAED) ,  Philadelphia,  over 
a  period  that  encompassed  the  27  to  29  March  1984  northeaster  (Gebert,  Watson, 
and  Hemsley,  in  preparation).  The  data  set  includes  profile  surveys  of  the 
subaerial  beach  taken  1  to  2  days  prior  to  arrival  of  the  storm  and  3  to 

4  days  after  the  storm.  Prestorm  subaqueous  profile  data  used  for  simulations 
were  collected  in  early  January  1984  and  poststorm  nearshore  surveys  were  made 
approximately  2  weeks  after  the  event.  A  time  series  of  wave  height  and 
period  was  recorded  by  a  wave  rider  buoy  located  directly  offshore  in  15. 2 -m 
water  depth,  and  water  level  at  Manasquan  Inlet  was  also  available. 

129.  Profile  measurements  were  made  along  eight  beach  transects  at 
Manasquan  to  the  north  of  the  Manasquan  Inlet  jetties  and  along  nine  transects 
at  Point  Pleasant  Beach  to  the  south  of  the  jetties  (Figure  10).  Point 
Pleasant  Beach  is  relatively  wide  and  sandy  with  no  coastal  structures  in  the 
surf  zone  south  of  the  678 -m- long  south  jetty  at  Manasquan  Inlet,  whereas 
Manasquan  has  a  narrow  sandy  beach  with  short  groins  located  at  approximately 
200-  to  300-m  intervals. 

130.  Regular  profile  surveys  were  carried  out  at  Manasquan  and  Point 
Pleasant  Beach  from  1982  to  1984  under  the  auspices  of  the  CE  Monitoring 
Completed  Coastal  Projects  (MCCP)  Program.  The  main  purpose  of  the  project  was 
to  monitor  the  structural  stability  of  the  jetties  after  rehabilitation  with 
16- ton  dolosse  armor  units.  Two  secondary  concerns  dealt  with  maintenance  of 
inlet  cross-sectional  area  and  potential  change  on  adjacent  beaches  caused  by 
the  jetties  (Gebert,  Watson,  and  Hemsley,  in  preparation).  Figure  10  gives  the 
locations  of  the  survey  lines  and  the  numbering  convention.  The  lines  are 
numbered  consecutively  from  1  starting  at  the  jetties  and  increasing  north  and 
south,  where  the  symbols  M  and  PP  are  used  to  denote  Manasquan  and  Point 
Pleasant  Beach,  respectively.  A  tide  gage,  operated  by  the  National  Oceano¬ 
graphic  Service  (NOS),  is  located  on  the  bay  side  of  Manasquan  Inlet.  At  the 
time,  a  wave  rider  buoy  was  operated  by  MCCP  and  located  in  15.2  m  of  water 
seaward  of  Manasquan.  Wave  conditions  generally  were  measured  every  6  hr  for  a 
20-min  period,  and  the  sampling  interval  was  automatically  decreased  to  1  or 

2  hr  when  the  wave  height  exceeded  3  m. 
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Figure  10.  Location  map  showing  profile  survey  lines  and  the  tide  gage 
and  wave  rider  buoy  positions  at  the  Manasquan  and  Point  Pleasant  Beach 

study  area 


131.  Toward  the  end  of  March  1984,  a  northeaster  passed  the  area 
causing  severe  erosion  of  the  beaches.  The  time-history  of  measured  signifi¬ 
cant  wave  height  and  peak  spectral  wave  period  for  the  storm  is  displayed  in 
Figure  11,  and  water  level  is  shown  in  Figure  12.  Water  surface  elevation  is 
given  with  respect  to  mean  sea  level  (MSL)  and  comprises  both  storm  surge  and 
tidal  variations.  Wave  height  began  to  increase  in  the  morning  on  28  March 
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and  peaked  about  24  hr  later  on  29  March  with  a  maximum  height  of  6.6  m.  By 
noon  on  30  March,  the  storm  had  passed,  and  more  typical  wave  heights  were 
recorded.  Wave  period  increased  sharply  with  arrival  of  the  storm,  similar  to 
the  behavior  of  wave  height,  but  long  periods  persisted  after  the  storm, 
creating  mild-steepness  waves  and  strongly  accretionary  conditions. 

132.  The  water  level  record  displayed  a  pronounced  peak  for  the  middle 
of  the  storm,  approximately  coinciding  with  the  maximum  wave  height.  The 
tidal  variation  at  Manasquan  Inlet  is  semidiurnal,  with  two  high  and  two  low 
waters  in  a  tidal  day,  and  the  mean  tidal  range  is  1.2  m.  Maximum  surface 
elevation  measured  during  the  storm  event  was  1.90  m  at  about  0600  on  29  March. 
Profile  measurements 

133.  Prestorm  beach  profile  surveys  were  taken  on  26  and  27  March  for 
Manasquan  and  Point  Pleasant  Beach,  respectively.  However,  when  the  storm 
arrived,  the  scheduled  Fathometer  surveys  had  to  be  halted  due  to  high  waves, 
and  only  wading  depth  profiles  were  obtained  for  these  dates.  Subaqueous 
profile  measurements  obtained  from  the  previous  profile  survey  (28  December 
1983  at  Manasquan  and  29  December  1983  at  Point  Pleasant  Beach)  were  used  to 
construct  complete  prestorm  profiles.  Poststorm  profiles  were  recorded  on  3 
April  at  Manasquan  and  2  April  at  Point  Pleasant  Beach.  Profile  survey 
elevations  were  made  with  respect  to  the  National  Geodetic  Vertical  Datum 
(NGVD) ,  and  were  converted  to  MSL  in  this  report  by  the  relationship  MSL(m)  - 
NGVD(m)  +  0.24  m. 

134.  Profile  survey  data  were  analyzed  to  yield  representative  profiles 
to  be  used  in  numerical  simulations.  Also,  longshore  variability  in  eroded 
volume  was  investigated  to  evaluate  the  applicability  of  a  two-dimensional 
approach.  Since  the  shoreline  was  curved  close  to  the  jetties,  especially  at 
Point  Pleasant  Beach,  profiles  at  the  respective  beaches  were  shifted  across  - 
shore  and  aligned  with  reference  to  MSL.  A  single  representative  profile  for 
each  beach  was  determined  by  averaging  profile  elevations  at  fixed  intervals. 

135.  Figure  13  shows  plots  of  the  eight  prestorm  profiles  measured  at 
Point  Pleasant  Beach.  The  shapes  of  the  profiles  are  surprisingly  similar, 
indicating  relatively  good  two -dimensionality.  However,  the  two  profiles 
located  closest  to  the  jetty  (PP-1  and  PP-2)  show  depths  that  are  signifi¬ 
cantly  less  than  those  associated  with  the  other  profiles.  This  shallower 
subaqueous  portion  of  the  profile  is  probably  a  result  of  impoundment  of  sand 
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Figure  13.  Profiles  measured  at  Point  Pleasant  Beach  during 
December  1983  aligned  with  respect  to  MSL 


transported  alongshore.  The  two  profiles  closest  to  the  jetty  were  therefore 
excluded  when  determining  the  generic  profile  for  simulation.  These  two 
profiles  also  exhibited  markedly  less  erosion  during  the  storm  as  compared 
with  the  other  profiles. 

136.  In  a  manner  similar  to  the  analysis  for  Point  Pleasant  Beach,  two 
profiles  were  removed  from  the  analysis  at  Manasquan;  seven  profiles  were  used 
for  determining  a  generic  profile  for  the  simulation.  The  transect  closest  to 
the  jetty  (M-l)  was  not  used  along  with  the  profile  farthest  to  the  north 
(M-9).  Their  exclusion  was  motivated  by  the  significantly  different  response 
of  the  two  profiles  to  the  storm,  for  which  substantial  accumulation  took 
place  at  their  seaward  ends.  This  is  attributed  to  longshore  effects  and  the 
influence  of  structures. 

137.  A  generic  prestorm  profile  was  created  using  the  average  profile 
from  the  immediate  prestorm  subaerial  surveys  matched  with  an  average 
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subaqueous  profile  developed  from  December  1983  survey  data.  In  the  same 
manner,  a  poststorm  profile  was  developed  based  on  the  surveys  of  2  and  3 
April  1984. 

138.  Representative  prestorm  and  poststorm  profiles  determined  for 
Manasquan  and  Point  Pleasant  Beach  (Figures  14  and  15)  differ  in  two  charac¬ 
teristics.  First,  the  dune  is  higher  and  narrower  at  Manasquan  as  compared 
with  Point  Pleasant  Beach.  Second,  a  wide  subaqueous  bar  exists  in  the 
profile  at  Manasquan,  extending  approximately  300  m  seaward  of  the  baseline, 
out  to  the  end  of  the  groins.  It  appears  as  if  the  groin  field  is  effective 
in  trapping  sand  on  the  nearshore  portion  of  the  profile. 

Beach  grain  size 

139.  Grain  size  is  a  principal  parameter  in  the  numerical  model  and  can 
be  difficult  to  quantify  for  use  in  SBEACH  since  a  wide  range  of  grain  sizes 
is  normally  present  on  a  natural  beach.  In  the  study  area,  investigations 
have  been  made  on  the  temporal  and  spatial  distribution  of  grain  size  (USAED, 
New  York  1954;  Gorman,  Stauble,  and  Reed,  in  preparation).  However,  because 
of  inconvenient  sample  location  and  lack  of  measurements  in  the  dune  region, 
available  grain  size  information  was  supplemented  with  sampling  from  the 
subaerial  portion  of  the  beach  for  specific  use  for  modeling  with  SBEACH. 

140.  The  additional  sediment  sampling,  performed  by  Mr.  Jeffrey  Gebert, 
USAED,  Philadelphia,  in  April  1989,  showed  a  uniform  median  grain  size  of 
approximately  0.5  mm  in  the  dune  region.  Median  grain  size  varied  more  on 
other  portions  of  the  subaerial  profile  with  coarser  material  typically 
present  at  the  berm  crest.  Sand  samples  were  taken  at  four  cross -shore 
locations  on  the  dry  portion  of  the  beach  and  at  four  alongshore  sites. 

Table  1  gives  a  summary  of  median  grain  size  data  available  for  the  Manasquan 
and  Point  Pleasant  Beach  study  sites. 

141.  Two  sampling  sites  south  of  the  jetties  and  one  to  the  north  were 
used  in  defining  median  grain  size  in  the  dune  region,  together  with  data  from 
previous  investigations.  Samples  collected  adjacent  to  the  jetties  were 
excluded  from  the  analysis  by  reasoning  similar  to  that  used  in  determining  a 
representative  profile.  An  average  median  grain  size  for  the  dune  region  was 
determined  to  be  0.49  mm  for  Manasquan  and  0.54  mm  for  Point  Pleasant  Beach. 

On  the  foreshore,  coarser  material  was  found  near  the  berm  crest  with  median 
grain  sizes  around  0.3  to  1.5  mm.  For  the  subaqueous  portion  of  the  profile, 
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. Prestorm  - Poststorm 

Figure  14.  Generic  prestorm  and  poststorm  profile  for  Manasquan 


Elevation  (m) 


. Prestorm  - Poststorm 

Figure  15.  Generic  prestorm  and  poststorm  profile  for  Point  Pleasant  Beach 
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Table  1 

Median  Grain  Size  (mm)  for  Manas quan  and  Point  Plus sant„ Beach 


Profile 

Landward  Limit 
of  Dune 

Dune  Crest 

Berm  Crest 

Mid- Tide 

Nearshore’ 

M-l* ** 

0.38 

0.43 

0.20 

M-2++ 

-- 

0.47 

0.56 

0.86 

M-9** 

-- 

0.44 

0.38 

0.16 

M-9+ 

-- 

0,51 

0.60 

1.31 

0.21 

PP-1** 

Point  Pleasant 

Beach 

0.28 

0.34 

0.21 

PP-1+ 

-- 

0.53 

-- 

0.36 

0.29 

PP.1++ 

0.54 

0.55 

1.51 

0.51 

-- 

PP-3++ 

0.47 

0.53 

1.13 

0.82 

-- 

PP-6++ 

0.48 

-- 

1.05 

0.68 

-- 

PP-8** 

-- 

0.53 

0.64 

0.19 

PP-8+ 

.  _ 

0.66 

0.41 

0.47 

0.18 

*  Average  D50  on  subaqueous  profile. 

**  June  1953. 

+  May  1987 . 

++  April  1989. 

median  grain  diameter  from  previously  collected  sand  samples  showed  an  overall 
trend  of  decreasing  grain  size  with  distance  offshore.  The  average  median 
grain  size  for  all  subaqueous  samples  showed  D5Q  -  0.19  mm  for  Manasquan  and 
D50  -  0.22  mm  for  Point  Pleasant  Beach. 

142.  An  indirect  method  of  obtaining  a  representative  grain  size  for 
the  subaqueous  portion  of  the  profile  is  to  approximate  the  beach  profile  with 
the  equation  proposed  by  Dean  (1977)  and  solve  for  the  shape  parameter  A  . 
From  the  shape  parameter,  it  is  possible  to  determine  the  corresponding  grain 
size  based  on  curves  presented  by  Moore  (1982)  and  Dean  (1987,  1990).  This 
equation  (h  -  Ax^/^)  was  least-square  fitted  to  the  prestorm  profile  from  the 
January  surveys  to  indirectly  determine  a  representative  grain  size  for  the 
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January  surveys  to  indirectly  determine  a  representative  grain  size  for  the 
subaqueous  portion  of  the  profile.  In  this  equation,  the  shape  parameter  A 
is  an  empirical  function  of  grain  size.  The  effective  median  grain  diameter 
calculated  for  profiles  M-2  to  M-8  was  0.47  mm,  and  for  profiles  PP-3  to  PP-8 
the  effective  grain  size  was  0.51  mm.  In  the  least-square  fit,  both  the  shape 
parameter  and  the  horizontal  location  of  the  still-water  shoreline  were 
estimated  to  minimize  the  sum  of  squares  between  measured  profiles  and  the 
profile  equation. 

143.  The  profile  equation  could  in  general  explain  a  large  portion  of 
the  variation  in  the  measurements,  but  visual  inspection  showed  some  deficien¬ 
cies  in  the  description.  To  illustrate  this  point,  Figure  16  shows  the  repre¬ 
sentative  prestorm  profile  at  Point  Pleasant  Beach  and  the  least-square  fitted 
profile.  For  an  average  median  grain  size  of  0.5  mm,  the  fitted  profile  gives 
an  overall  reasonable  description  relative  to  the  representative  profile. 
Directly  seaward  of  the  shoreline,  however,  the  representative  profile  is 
steeper  than  the  fitted  profile,  whereas  the  opposite  occurs  in  the  offshore 
region.  The  actual  profile  shape  is  probably  sensitive  to  varying  median 
grain  size,  resulting  in  a  shape  that  locally  deviates  from  a  fit  involving  a 
2/3  power  for  a  single  grain  size. 

Alangshpye  variability  in  gtorsLergsion 

144.  The  numerical  model  simulates  cross -shore  sand  transport  only; 
longshore  effects  are  not  represented.  In  order  to  evaluate  this  approxi¬ 
mation  for  the  target  storm  event,  the  amount  of  erosion  was  determined  for 
each  profile  survey.  Equal  amounts  of  erosion  would  suggest  that  the  forcing 
conditions  were  similar  during  the  storm  event  and  differentials  in  longshore 
transport  were  small.  It  is  important  to  note  that  significant  recovery 
occurred  before  the  poststorm  profile  was  surveyed.  Figure  17  illustrates  the 
recorded  volume  difference  above  the  1-m  contour  between  the  prestorm  and 
poststorm  profile.  Profile  location  is  referenced  to  the  profile  closest  to 
the  jetties,  going  north  and  south  at  Manasquan  and  Point  Pleasant  Beach, 
respectively.  Table  2  lists  measured  volumes  of  erosion  above  the  1  and  0 
(MSL  shoreline)  contours. 

145.  Significantly  less  erosion  is  indicated  for  the  two  profile  lines 
closest  to  the  jetties  at  Point  Pleasant  Beach.  South  of  profile  line  PP-2, 
the  eroded  volume  is  fairly  constant  (50  to  65  m^/m)  except  for  line  PP-6. 
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. Prestorm  - Least-square 

Figure  16.  Generic  prestora  profile  at  Point  Pleasant  Beach 
and  least-square  fitted  equilibrium  profile 


Point  Pleasant  Beach  Manasquan 

Figure  17.  Volume  difference  above  the  1-m  contour  between  prestorm 
and  poststorm  profiles  at  Manasquan  and  Point  Pleasant  Beach 
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Table  2 

Eroded  Profile  Volumes  at  Manasauan  and  Point  Pleasant  Beach 
for  the  March  1984  Storm  Event  (m3/!!!) 


Contour 

Grouo  Averaee 

Average  of  Selected 
Profiles 

Minimum 

Maximum 

0 

36.3 

Manasauan 

39.3 

14.7 

62.6 

1-m 

32.4 

34.7 

13.4 

56.2 

0 

Point  Pleasant  Beach 

48.9  55.8 

9.6 

66.0 

1-m 

45.8 

53.8 

8.5 

56.2 

Examination  of  the  poststorm  profile  survey  at  line  PP-6  shows  considerable 
deposition  on  top  of  the  dune.  As  previously  noted,  the  eroded  volume 
calculation  interval  includes  a  period  of  recovery  or  accretion  that  contami¬ 
nates  the  data  set  with  respect  to  storm- induced  erosion  modeling.  The 
average  volume  eroded  above  the  1-m  contour  for  all  profile  lines  at  Point 
Pleasant  Beach  is  45.8  m^/m,  and  for  profiles  PP-3  to  PP-8  the  average  eroded 
volume  is  53.8  m^/m.  To  gage  the  intensity  of  storm  activity  and  resultant 
magnitude  of  beach  profile  change,  data  from  this  event  were  compared  with 
similar  parameters  from  13  other  east  coast  storms  (Birkemeier,  Savage,  and 
Leffler  1988).  Comparison  of  peak  wave  and  water  level  conditions,  and 
associated  changes  in  beach  profile  volume,  with  data  from  the  27  to  29  March 
1984  event  at  Manasquan  and  Point  Pleasant  Beach  indicates  that  dynamic  charac¬ 
teristics  associated  with  this  storm  and  resultant  volume  changes  were  more 
extreme  than  any  of  the  13  previously  recorded  events.  Analysis  of  extreme 
wave  statistics  for  this  area  suggests  that  this  major  storm  event  has  a 
frequency  of  occurrence  of  about  30  years  (Jensen  1983,  Corson  and  Tracy  1985). 

146.  The  profiles  at  Manasquan  experienced  less  erosion  than  those  at 
Point  Pleasant,  with  an  average  eroded  volume  for  all  profile  lines  of 
32.4  m^/m.  Survey  line  M-l  closest  to  the  jetty  experienced  somewhat  less 
erosion  than  all  but  one  of  the  other  profiles.  The  average  volume  eroded 
above  the  1-m  contour  for  profiles  M-2  to  M-8  was  34.7  m^/m. 
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147.  Variability  in  volume  of  erosion  is  greater  for  Manasquan  than  for 
Point  Pleasant  if  the  profile  lines  closest  to  the  jetty  are  neglected 
(Figure  17) .  This  is  probably  due  to  localized  processes  associated  with  the 
groin  field,  giving  rise  to  longshore  variability  in  sand  transport.  The 
lower  average  eroded  volume  for  Manasquan  is  also  partially  attributed  to  the 
groins,  which  evidently  caused  a  large  subaqueous  bar  to  develop.  This 
feature  increases  energy  dissipation  on  the  seaward  portion  of  the  profile  and 
reduces  wave  energy  reaching  the  dune  face.  Expressed  differently,  the  bar  is 
a  sand  reservoir  that  must  first  be  depleted  by  erosive  waves  before  they  can 
fully  attack  the  subaerial  beach. 

148.  In  summary,  it  appears  that  longshore  variability  in  eroded  volume 
was  relatively  small  for  Point  Pleasant  Beach  because  structures  were  not 
present  to  create  differentials  in  longshore  transport.  At  Manasquan,  larger 
volume  differences  are  noted  due  to  the  presence  of  structures  and  the  location 
of  profile  lines  with  respect  to  the  groins.  Applying  the  numerical  model  to  a 
representative  generic  profile  should  be  acceptable  both  at  Manasquan  and  Point 
Pleasant  Beach  in  terms  of  longshore  variability  in  eroded  volume.  However,  if 
only  one  profile  survey  were  available  for  modeling,  significant  variability  or 
deviation  in  calculated  and  measured  results  could  be  expected. 

Numerical  simulations 

149.  Average  representative  profiles  were  used  in  the  numerical  simula¬ 
tions.  Since  median  grain  diameter  did  not  differ  significantly  between  the 
two  sites,  the  same  value  was  specified.  A  schematized  description  of  cross¬ 
shore  variation  in  grain  size  was  specified  based  on  sampling  results,  with  a 
diameter  of  0.5  mm  shoreward  of  the  MSL  shoreline  and  0.3  mm  seaward  of  the 
MSL  shoreline.  Model  parameters  were  initially  set  at  the  same  values  as  for 
the  large  wave  tank  studies,  but  some  modifications  were  made  in  the  calibra¬ 
tion  process. 

150.  Simulation  periods  for  the  two  beaches  corresponded  with  different 
survey  times  at  the  respective  sites.  The  exact  time  of  day  of  the  surveys 
was  not  known  and  was  assumed  to  be  1200.  For  Manasquan  the  simulation 
interval  was  26  March  to  3  April,  and  for  Point  Pleasant  Beach  it  was  27  March 
to  2  April.  Wave  and  water  level  conditions  for  the  calculation  periods  are 
given  in  Figures  11  and  12,  for  which  cubic  spline  interpolation  was  used  to 
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obtain  intermediate  values  between  measurement  points.  The  model  time  step 
was  5  min,  and  the  length  step  was  2  m. 

151.  Simulation  results  are  shown  in  Figure  18  for  Manasquan  and  in 
Figure  19  for  Point  Pleasant  Beach.  The  initial  and  poststorm  representative 
profiles  are  plotted  with  the  calculated  profiles.  Calculated  erosion  was 
defined  as  the  time  at  which  the  model  switched  from  calculating  erosion  to 
accretion  by  Equation  1,  and  recovery  was  defined  as  the  time  at  which  the 
poststorm  survey  was  taken.  Since  poststorm  profile  measurements  were  made 
after  several  days  of  recovery  waves,  the  final  measured  profile  includes  a 
period  of  potentially  significant  beach  recovery.  In  particular,  at  Point 
Pleasant  the  recovery  was  very  strong  with  a  wide  gentle  ridge  evolving  with  a 
crest  2  m  above  the  MSL  shoreline.  Beach  recovery  at  Manasquan  was  consider¬ 
able  but  not  as  pronounced  as  at  Point  Pleasant  Beach,  probably  due  in  part  to 
a  steeper  dune  face  and  a  narrower  foreshore,  conditions  that  tend  to  retard 
the  recovery  process. 

152.  Manasquan.  Simulated  profile  retreat  at  the  top  of  the  dune  is 
close  to  the  measured  retreat  for  Manasquan,  and  the  dune  slope  is  reproduced 
well.  Also,  the  poststorm  eroded  volume  is  considered  to  be  in  good  agreement 
with  what  probably  occurred  as  inferred  from  the  profile  taken  on  3  April  by 
extrapolating  the  slope  of  the  dune  face  to  MSL  and  neglecting  the  accretionary 
berm  feature.  Model  simulation  of  recovery  is  limited  to  qualitative  and  not 
quantitative  agreement  with  the  measurements,  showing  berm  buildup  but  with  a 
smaller  volume.  Another  aspect  of  profile  development  that  differs  between 
calculated  and  measured  is  the  evolution  of  the  bar  located  near  the  seaward 
terminus  of  the  groins.  The  numerical  model  predicts  that  the  upper  part  of 
the  terrace  is  transported  farther  offshore,  whereas  measurements  indicated 
that  this  feature  was  stable  during  the  event.  This  shortcoming  of  the  model 
may  be  partially  attributed  to  neglect  of  longshore  processes. 

153.  Point  Pleasant  Beach.  The  model  simulation  for  Point  Pleasant 
Beach  shows  much  the  same  characteristics  as  for  Manasquan,  although  predicted 
dune  retreat  was  smaller  as  compared  with  measured  retreat.  Poststorm  eroded 
volume  is  again  considered  to  be  in  good  agreement  with  the  actual  amount 
through  examination  of  the  2  April  profile  and  disregarding  the  large  ridge 
assumed  to  have  beei.  formed  during  the  recovery  phase.  Simulated  recovery  is 
stronger  at  Point  Pleasant  compared  with  Manasquan,  which  is  qualitatively  in 
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Figure  18.  Measured  and  calculated  profile  change  at  Manasquan 
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Figure  19.  Measured  and  calculated  profile  change  at 
Point  Pleasant  Beach 
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accordance  with  observations,  although  the  predicted  berm  volume  is  smaller. 
The  numerical  model  also  indicates  development  of  bar-like  features  in  the 
seaward  portion  of  the  profile  that  were  not  present  in  the  measured  profile. 
Discussion 

154.  In  the  simulations  for  both  beaches,  the  transport  rate  coeffi¬ 
cient  K  was  changed  to  obtain  optimal  agreement  between  model  results  and 
the  measured  profile.  The  coefficient  value  determined  for  both  sites  was 
K  -  2.0  x  10 m^/N  ,  compared  with  the  value  from  the  large  wave  tank 
simulations  of  K  -  1.6  x  10*^  m^/N  and  from  the  Duck,  NC,  simulations  of 

K  -  0.7  x  10 m^/N  .  The  Manasquan  and  Point  Pleasant  Beach  profile  change 
simulations  used  an  optimal  K- value  similar  to  the  large  wave  tank  simulation, 
whereas  the  coefficient  for  the  Duck  data  was  considerably  lower,  probably  due 
to  the  unusual  stability  of  the  shoreline  at  Duck  that  hindered  erosion. 

155.  To  better  account  for  naturally  variable  wave  conditions,  wave 
heights  used  in  model  simulations  were  randomly  varied  a  certain  percentage 
around  the  values  interpolated  from  the  measured  wave  heights.  This  technique 
simulates  variability  in  the  location  of  the  break  point  and  development  of 
profile  features.  The  method  for  randomly  varying  wave  height  consisted  of 
choosing  a  height  from  a  rectangular  distribution  with  its  center  as  the 
interpolated  height  and  its  width  as  a  predefined  percentage  of  the  interpo¬ 
lated  height.  Different  distribution  widths  were  tested,  and  10  to  20  percent 
proved  to  give  a  good  qualitative  description  of  profile  evolution.  For  the 
calibration  runs  displayed  in  Figures  17  and  18,  a  20-percent  variation  was 
employed.  Although  this  technique  was  not  performed  in  simulations  for  the 
case  at  Duck,  a  larger  K- value  and  smoother,  more  realistic  profile  features 
would  be  expected. 

156.  In  order  to  obtain  a  qualitatively  improved  description  of  the 
accretion  process,  adjustments  to  the  transport  rate  distribution  were  made. 
For  the  large  wave  tank  data  simulations,  the  transport  rate  distribution  was 
calculated  the  same  way,  independent  of  whether  onshore  or  offshore  transport 
occurred.  However,  for  a  natural  beach  with  .■?  variable  wave  climate,  unusual 
conditions  arise  that  must  be  dealt  with  in  a  special  way.  For  example,  some 
acc^etionary-type  waves  are  small  enough  that  breaking  will  not  occur  until 
they  arrive  almost  at  the  shoreline.  The  transport  rate  distribution  in  this 
case  is  very  narrow  and  peaked,  producing  large  accretion  in  the  cells  around 
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the  shoreline.  This  local  accumulation  steepens  the  lower  foreshore  slope, 
making  more  waves  break  at  the  shoreline.  To  take  care  of  this  condition,  if 
at  any  time-step  the  surf  zone  was  calculated  to  be  less  than  four  cells  wide, 
the  transport  rate  was  specified  as  linearly  decreasing  from  the  plunge  point 
to  the  runup  limit.  The  transport  rate  at  the  plunge  point  was  determined 
from  the  wave  energy  dissipation  according  to  Equation  4,  as  usual. 

157.  If  onshore  transport  is  high,  berm  buildup  is  rapid,  and  the 
foreshore  beach  slope  may  become  unrealistically  steep,  limited  only  by  the 
angle  of  initial  yield.  Such  slopes  predicted  by  the  numerical  model  are  in 
contradiction  with  observations  from  the  large  wave  tank  data  (Takeda  and 
Sunamura  1983,  Larson  and  Kraus  1989a)  and  much  steeper  than  what  is  found  on 
a  natural  beach  (Bascom  1951,  Wiegel  1964,  Thomas  and  Baba  1986,  Hughes  and 
Cowell  1987).  To  limit  steepening  of  the  foreshore,  the  cross-shore  sand 
transport  rate  was  reduced  with  respect  to  the  local  slope  when  transport  was 
directed  onshore.  This  introduced  a  second  dependence  on  beach  slope  because 
Equation  4  already  contains  a  term  dependent  on  dh/dx  ,  but  was  necessary  to 
achieve  a  smoother,  more  realistic  approach  to  equilibrium. 

158.  Eroded  volumes  for  the  predicted  poststorm  profile  and  profile 
after  recovery  were  calculated  above  the  0  and  1-m  contour  and  are  displayed 
In  Table  3.  These  values  can  be  compared  with  eroded  volumes  based  on 
measured  profile  evolution  listed  in  Table  2.  As  previously  discussed,  eroded 


Table  3 

Calculated  Eroded  Profile  Volumes  at  Manasauan 
and  Point  Pleasant  Beach  (m3/m^ 


Contour 

Maximum  Erosion 

With  Recovery 

Manas quan 

0 

51.5 

44.2 

1-m 

37.2 

37.2 

Point  Pleasant  Beach 

0 

86.4 

73.9 

1-m 

54.6 

54.2 
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volume  determined  from  the  measured  profiles  Includes  a  period  of  recovery  as 
the  Immediate  poststorm  profile  was  not  surveyed.  Since  model  prediction  of 
recovery  was  only  in  qualitative  agreement  with  measurements,  calculated 
eroded  volumes  at  the  end  of  the  simulation  period  were  too  large  in  com¬ 
parison  with  measured  changes.  Calculated  eroded  volume  at  Manasquan  agrees 
well  with  measurements,  although  the  shape  of  the  profile  differs,  with  the 
measured  profile  exhibiting  less  foreshore  erosion  and  greater  dune  retreat. 
Poststorm  erosion  is  probably  well  predicted  at  Point  Pleasant  Beach,  but  berm 
buildup  is  underestimated.  Little  or  no  recovery  is  produced  by  the  numerical 
model  above  the  1-m  contour,  which  differs  from  what  occurred  in  the  field. 

159.  In  conclusion,  the  numerical  model  was  successful  in  reproducing 
the  erosional  portion  of  the  studied  storm  event,  but  predicted  profile  evolu¬ 
tion  during  the  recovery  stage  was  only  in  qualitative  agreement  with  measure¬ 
ments.  Parameter  values  obtained  for  model  calibration  were  similar  to  those 
found  for  simulations  of  the  large  wave  tank  data,  particularly  the  transport 
rate  coefficient  K  .  A  two-dimensional  approach  to  estimate  profile  evolu¬ 
tion  during  a  storm  event  was  justified  based  on  the  small  measured  alongshore 
variability  in  eroded  volume.  However,  if  nearshore  structures  such  as  groins 
exist  in  a  study  area,  care  should  be  taken  in  choosing  a  representative 
profile,  and,  if  possible,  averages  of  several  profiles  should  be  used. 

Storffi.Evqnt  <?n  VS  w<»sp 

160.  A  mission-support  project  performed  by  CERC  for  the  Coast  of 
California  Storm  and  Tidal  Waves  Study  conducted  by  USAED,  Los  Angeles, 
provided  an  opportunity  to  apply  SBEACH  to  storm-erosion  events  on  the 
southern  coast  of  California  (Kadib,  Fine,  and  Ryan  1989).  This  application 
involved  considerably  different  environmental  conditions  than  east  coast  field 
tests  of  the  model.  West  coast  conditions  of  note  include  small  storm  surge, 
typically  less  than  0.3  m,  longer  period  waves  than  on  the  east  coast,  and 
broad,  low- relief  bars  if  bars  are  present. 

161.  Documentation  of  storm- induced  beach  erosion  suitable  for  model 
testing  and  calibration  is  relatively  scarce  for  the  west  coast,  and  no 
publications  relating  to  numerical  modeling  of  dune  erosion  were  found  in  a 
literature  review.  As  of  late  1988,  after  a  search  of  the  literature  and 
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inquiry  to  leading  coastal  investigators,  only  the  data  collected  by  Aubrey 
(1978)  were  considered  to  be  adequate.  The  measurements  of  interest  were  made 
at  Torrey  Pines  Beach,  California,  located  approximately  2  miles  north  of 
Scripps  Institution  of  Oceanography  (Aubrey,  Inman,  and  Nordstrom  1976;  Aubrey 
1978;  Gable  1979),  between  December  1975  and  January  1977. 

162.  Figure  20  shows  the  locations  of  three  range  lines  along  Torrey 
Pines  Beach  referred  to  as  North  Range,  Indian  Canyon  Range,  and  South  Range. 
Onshore  and  offshore  portions  of  the  profiles  were  measured  periodically 
between  1972  and  1978,  together  with  the  associated  wave  climate  and  water 
level.  All  elevations  were  measured  relative  to  NGVD;  MSL  is  4.3  cm  above 
NGVD  in  this  region.  Time  series  of  wave  height,  wave  period,  and  water  level 
for  the  periods  29  December  1975  to  9  January  1976  and  21  December  1976  to 

4  January  1977  were  available  for  defining  short-term  beach  profile  change 
produced  by  storms.  Wave  measurements  were  made  along  Indian  Canyon  Range 
with  an  array  of  pressure  sensors  located  at  a  depth  of  approximately  9.1m 
(Figure  20) .  Examination  of  the  Torrey  Pines  data  shows  these  events  to  be 
less  significant  than  most  storms  producing  notable  beach  change  in  southern 
California.  However,  this  information  provides  the  only  comprehensive  data 
base  on  storm  conditions  and  associated  beach  profile  response  for  the  region 
that  could  be  found  during  this  study. 

163.  Beach  profile  data  were  collected  three  times  during  the  first 
storm  event  (12/29/75,  12/31/75,  and  1/9/76)  and  twice  during  the  second  event 
(12/21/76  and  1/4/77).  All  surveys  were  carried  out  near  the  time  of  predicted 
low  water.  The  offshore  bathymetry  indicated  very  little  longshore  variation, 
implying  only  minor  differences  in  wave  properties  at  the  three  range  lines. 
Therefore,  an  average  profile  shape  was  constructed  for  each  date  by  matching 
the  zero  crossings  (referenced  to  NGVD)  for  each  of  the  three  profiles  and 
averaging  the  elevations  at  known  positions  relative  to  z  -  0  .  Use  of  an 
average  profile  should  reduce  changes  associated  with  longshore  variability 
caused  by  local  differences  in  nearshore  bottom  topography,  sediment  size,  and 
waves  (Savage  and  Birkemeier  1987,  Birkemeier  et  al.  1988).  These  composite 
data  were  used  for  simulating  beach  profile  response  to  both  storm  events.  The 
recorded  time  series  of  waves  and  water  level  between  surveys  were  input  to  the 
beach  profile  change  model  SBEACH,  and  data  quality  checks  were  made  on  all 
process  and  response  measurements  before  model  runs  proceeded. 
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Figure  20.  Location  map  of  Torrey  Pines  Beach  showing  position 
of  range  lines  and  pressure  sensor  array  (after  Aubrey  1978) 


164.  Model  parameter  values  resulting  from  calibration  with  the  large- 
scale  laboratory  data  sets  (Report  1)  were  initially  used  to  simulate  field 
beach  profile  response.  It  became  apparent,  however,  that  values  of  the  two 
principal  empirical  coefficients  would  have  to  be  modified  to  achieve  optimal 
agreement  between  measured  and  calculated  profiles.  The  sand  transport  rate 
coefficient  K  in  Equation  4  obtained  from  these  two  west  coast  storm  events 
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was  slightly  smaller  than  that  obtained  from  the  laboratory  experiments  (1.4  x 
10*^  m^/N  versus  1.6  x  10*^  m^/N) ,  and  the  transport  rate  coefficient  e  for 
the  slope -dependent  term  of  the  transport  rate  equation  (Equation  37)  was 
increased  from  0.001  to  0.002.  In  addition,  to  account  for  time-varying, 
random  wave  height  fluctuations  that  tend  to  create  relatively  smooth  macro¬ 
scale  features  observed  on  the  subaqueous  profile  in  nature,  a  30-percent 
random  variation  in  recorded  wave  height,  representing  one  standard  deviation 
of  the  measured  wave  height  for  each  6 -hr  time  interval,  was  imposed  on  the 
entire  time  series.  A  time  step  of  5  min  and  a  length  step  of  2  m  were  used 
for  all  simulations. 

165.  Finally,  it  was  inferred  that  water  level  measured  at  the  gage 
(located  approximately  1  km  offshore)  was  not  entirely  representative  of  water 
elevations  at  the  shoreline  influencing  profile  change.  Calculated  beach 
profile  response  best  matched  measured  profile  change  when  water  level  at  the 
gage  was  increased  by  approximately  0.3  m  during  the  storm  event.  This  is 
consistent  with  onshore-  and  offshore-directed  wind  stress  associated  with 
storms  affecting  the  southern  California  coast  that  typically  create  fluctua¬ 
tions  in  water  surface  elevation  at  the  shoreline,  producing  setup  as  a  storm 
approaches  and  setdown  as  a  storm  passes. 

Storm  Event  1 

166.  Figure  21  shows  the  time  series  of  wave  height,  peak  wave  period, 
and  water  level  for  the  period  29  December  1975  through  9  January  1976.  These 
measured  quantities  were  available  at  6-hr  intervals.  Pressure  data  from  the 
array  of  sensors  were  depth-corrected  and  converted  to  sea-surface  elevation, 
yielding  the  significant  wave  height.  Since  the  time  step  used  in  the 
numerical  model  (typically  5  min)  was  shorter  than  the  wave  and  water  level 
measurement  interval,  input  values  for  time  steps  between  measurements  were 
interpolated. 

167.  A  2-day  period  near  the  beginning  of  the  record  illustrates  storm 
energy  resulting  from  a  relatively  weak  North  Pacific  low  pressure  system. 
Maximum  wave  heights  «ere  approximately  1.8  m  with  a  peak  period  of  11.1  sec. 
Elevated  water  levels  associated  with  storm  surge,  coincident  with  high  tide, 
produced  a  measured  water  level  at  the  gage  of  1.2  m  at  the  time  of  maximum 
wave  height.  Poststorm  wave  height  decreased  to  less  than  0.9  m,  and  peak 
wave  period  increased  as  swell  wave  energy  became  dominant. 
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168.  The  distribution  of  median  grain  diameter  along  the  profile  is 
required  to  run  SBEACH.  Grain  size  data  were  not  available  for  this  period, 
but  median  grain  diameter  recorded  near  the  site  in  October  1983  showed  little 
variation  between  the  1-  and  -1-m  contours  (USAED,  Los  Angeles  1984). 

Therefore,  an  average  median  grain  diameter  of  0.17  mm  was  used  for  model 
calibration. 

169.  Since  this  storm  was  of  minimal  intensity,  recorded  beach  profile 
erosion  was  relatively  minor.  Figure  22  illustrates  measured  and  calculated 
profile  shape  for  the  11 -day  period.  Calculated  profile  change  from 

29  December  to  31  December  1975  was  concentrated  on  the  upper  half  of  the 
foreshore  and  was  related  to  locally  elevated  water  levels.  Generally, 
calculated  profile  response  above  NGVD  was  slightly  less  than  measured 
adjustments;  the  2-m  contour  retreated  approximately  4.6  m,  and  profile  volume 
change  above  the  mean  shoreline  (z  -  0)  was  -8  m^/m.  Only  small  adjustments 
in  bar  elevation  were  recorded. 

170.  Measured  profile  response  between  29  December  1975  and  9  January 
1976  includes  effects  of  the  2-day  storm  event  as  well  as  approximately  8  days 
of  poststorm  recovery  wave  action.  After  31  December  1975,  most  profile 
changes  occurred  below  NGVD,  where  two  new  bar  features  developed.  However, 
cumulative  retreat  of  the  2-m  contour  increased  to  10  m,  and  associated 
subaerial  volume  change  increased  to  -15  m^/m.  The  calculated  profile  for 

9  January  1976  shows  a  small  bar  located  at  about  the  0.6-m  depth  contour,  in 
close  proximity  to  the  innermost  measured  bar.  Although  the  measured  bar 
immediately  seaward  of  this  position  is  not  represented  on  the  calculated 
profile,  the  broad,  low-relief  deposit  at  about  the  -2.4-m  contour  is  fairly 
well  represented  by  the  simulation.  Calculated  profile  change  above  NGVD 
showed  increased  retreat  of  the  2-m  contour  between  31  December  1975  and 
9  January  1976,  resulting  in  a  cumulative  retreat  of  about  6.4  m.  Simulated 
profile  volume  adjustments  gave  a  cumulative  loss  of  14.8  m^/m.  A  comparison 
of  measured  and  calculated  profile  changes  shows  good  agreement  for  this  storm. 
Storm  Event  2 

171.  Verification  of  SBEACH  simulations  was  accomplished  using  process 
and  response  data  from  a  second  small  storm  event  that  occurred  at  Torrey 
Pines  Beach  approximately  1  year  after  Storm  Event  1.  Data  were  analyzed  and 
prepared  for  model  application  in  a  manner  identical  to  that  for  the  winter 
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Figure  22.  Measured  and  calculated  beach  profile  change, 
Torrey  Pines,  winter  1975 


1975  storm.  Figure  23  shows  the  time  series  of  wave  and  water  level  para¬ 
meters  for  the  period  21  December  1976  to  4  January  1977.  In  this  case,  the 
record  covers  14  days  with  an  increase  in  wave  height  occurring  toward  the 
middle  portion  of  the  time  series.  Small  variations  in  wave  period  suggest 
that  locally  generated  winds  likely  were  not  the  main  source  of  energy  for 
this  event.  Consequently,  an  increase  in  water  level  due  to  wind  setup  at  the 
shoreline  was  not  imposed,  and  only  direct  measurements  of  water  level  at  the 
gage  were  used  for  simulations.  All  other  previously  determined  model 
parameters  were  held  constant. 

172.  A  maximum  wave  height  of  1.6  m  was  recorded  on  27  December  1976 
associated  with  a  15.1 -sec  wave  period.  For  the  next  4.5  days,  wave  height 
averaged  about  1.1  m.  Fluctuations  in  water  level  were  at  a  minimum  during 
this  time;  however,  water  surface  elevation  remained  above  NGVD  while  wave 
heights  were  greatest.  This  feature  was  not  present  in  Storm  Event  1  and 
likely  contributed  to  an  increase  in  eroded  volume  above  NGVD  during  Storm 
Event  2.  Wave  energy  associated  with  North  Pacific  swell  was  dominant  in  the 
14 -day  record. 

173.  This  second  event  exhibited  lower  wave  energy  than  the  first  storm; 
however,  changes  in  profile  volume  and  contour  position  above  NGVD  were 
slightly  greater.  This  can  be  attributed  to  the  duration  of  water  level  above 
the  mean  water  line  during  higher  wave  conditions.  Figure  24  illustrates 
measured  and  calculated  profile  change  for  the  period  21  December  1976  to 

4  January  1977.  Unlike  responses  recorded  during  the  first  storm  event,  a 
relatively  thin  veneer  of  sand  was  eroded  from  the  original  profile  between  the 
berm  crest  (approximately  2.3-ra  contour)  and  the  -2.1-m  contour.  This  volume 
of  sand  was  not  conserved  along  the  profile  since  deposits  of  near-equal  volume 
were  not  recorded  at  a  seaward  position.  Redistribution  of  sand  by  longshore 
processes  is  a  likely  cause  for  the  lack  of  cross -shore  conservation  of  mass 
during  this  period  of  dominant  swell  energy.  Measurements  of  beach  profile 
change  above  NGVD  showed  a  maximum  retreat  associated  with  the  0-depth  contour 
(12.8-m  recession)  and  a  steady  decrease  in  contour  retreat  at  higher  eleva¬ 
tions.  The  1-m  contour  receded  8.8  m,  and  the  2-m  contour  receded  4  ra. 

Measured  changes  in  profile  volume  above  NGVD  were  approximately  -16.3  m^/m. 

174.  Calculated  profile  change  satisfactorily  verified  predictive 
capabilities  of  the  model.  Changes  in  position  of  the  1-  and  2-m  contours 
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Figure  23.  Time  series  of  wave  and  water  level  characteristics, 

winter  1976  storm 
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Figure  24.  Measured  and  calculated  beach  profile  change, 

Torrey  Pines,  winter  1976 

were  nearly  identical  to  measured  changes  (-10.4  and  -4  m) ,  and  the  0-depth 
contour  was  slightly  underpredicted.  Therefore,  profile  volume  change  above 
NGVD  was  consistent  with  measured  changes  (-17.6  mV®) .  However,  calculated 
profile  shape  on  the  seawardmost  portion  of  the  active  profile  was  signifi¬ 
cantly  different  from  measured  change.  This  apparent  inconsistency  is 
explained  by  examining  the  way  in  which  profile  change  is  calculated.  At  the 
shoreward  end  of  the  profile,  the  runup  limit  constitutes  a  boundary  where 
cross-shore  sand  transport  becomes  zero.  The  seaward  boundary  is  determined 
by  the  depth  at  which  no  significant  sediment  motion  occurs.  Once  the  distri 
bution  of  the  sand  transport  rate  is  known,  profile  change  is  calculated  from 
the  mass  conservation  equation.  If  the  assumption  of  conservation  of  mass  in 
a  two-dimensional  sense  is  not  approximated  by  a  three-dimensional  natural 
process,  either  an  excess  or  deficit  of  sand  volume  will  be  predicted.  In 
this  case,  calculated  profile  change  shows  an  excess  of  sand  volume  deposited 
as  a  bar,  suggesting  that  longshore  processes  may  have  been  important  in 
redistributing  nearshore  sand  deposits.  Verification  of  calculated  profile 
response  for  this  region  was  otherwise  good. 
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PART  VI:  MODEL  APPLICATION  TO  BEACH-FILL  DESIGN 


175.  Beach  nourishment  is  becoming  a  preferred  shore  protection  measure 
because  it  is  often  the  least  expensive  solution,  is  flexible  and  can  augment 
existing  shore  protection  structures,  and  provides  a  natural  and  enjoyable 
coastal  environment.  Prediction  of  beach- fill  behavior  is  required  to 
determine  the  benefits  of  a  nourishment  project  in  relation  to  other  shore 
protection  alternatives  as  well  as  to  estimate  its  longevity  and  function¬ 
ality.  A  realistic  description  of  coastal  sediment  transport  processes  and 
beach- fill  change  involves  numerous  hydrodynamic,  sediment,  and  morphologic 
characteristics  that  vary  in  space  and  time.  Numerical  modeling  of  beach 
profile  and  fill  change  provides  a  general  and  quantitative  means  of  repre¬ 
senting  the  combined  and  simultaneous  effects  of  these  variables.  Application 
to  beach- fill  design  was  one  of  the  objectives  in  developing  SBEACH. 

176.  Figure  25  is  a  definition  sketch  of  a  nearshore  region  containing 
a  beach  fill.  The  fill  is  acted  on  by  ordinary  cross-shore  wave  processes, 
such  as  those  associated  with  episodic  and  potentially  catastrophic  storms, 
which  mold  it  from  an  initial  shape  to  a  dynamic  equilibrium  shape.  The  fill 


Figure  25.  Definition  sketch  categorizing  a  beach  fill 
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is  also  exposed  to  longshore  currents  induced  by  waves  and  wind,  which  move 
sand  alongshore  and  deplete  the  fill  when  updrift  sand  source  is  inadequate  or 
maintenance  nourishment  is  not  provided.  Longshore  transport  does  not  enter 
in  SBEACH;  therefore,  the  model  is  most  appropriately  applied  to  the  internal 
region  of  a  long  fill  and  far  from  coastal  structures  and  boundaries  where 
differentials  in  longshore  transport  may  be  large. 

177.  This  chapter  describes  hypothetical  applications  of  SBEACH  to 
predict  the  initial  adjustment  of  beach  fill  to  wave  action  by  storms. 
Initially,  simulations  are  made  for  two  types  of  storms  occurring  on  the  US 
east  coast.  One  synthetic  tropical  storm  (hurricane)  and  one  synthetic 
extratropical  storm  (northeaster)  representing  mid-Atlantic  events  with 
approximate  2-  to  5-year  return  periods  are  used  to  predict  erosion  associated 
with  two  hypothetical  beach-fill  configurations  and  subsequent  poststorm 
recovery.  Eroded  sand  volume  and  contour  movement  are  evaluated  as  a  function 
of  storm  type,  fill  cross  section,  grain  size,  and  time, 

178.  In  the  second  half  of  the  chapter,  simulations  are  made  of  the 
response  of  two  different  fill  cross  sections  of  two  different  grain  sizes  to 
a  tropical  storm  with  a  return  period  of  about  10  years.  Comparisons  are 
made  to  a  beach  without  a  fill  to  evaluate  the  effectiveness  of  artificial 
nourishment  as  a  means  of  shore  protection.  A  seawall  is  included  in  the 
simulations  to  illustrate  model  capabilities  to  predict  profile  evolution  in 
situations  involving  a  vert '.cal  wall.  Poststorm  recovery  is  also  simulated 
for  this  case  making  it  feasible  to  assess  the  possibility  for  the  eroded  fill 
sand  to  contribute  material  in  the  beach  recovery  process.  Some  general 
guidelines  are  provided  concerning  the  design  of  beach  nourishment  schemes 
based  on  the  numerical  simulations. 

179.  Finally,  a  concluding  discussion  is  given  regarding  capabilities 
and  limitations  of  the  numerical  model  as  a  tool  for  predicting  beach- fill 
evolution.  Much  of  the  material  presented  in  this  chapter  may  be  found  in 
Kraus  and  Larson  (1989)  and  Larson  and  Kraus  (in  press). 

Fill  Response  to  Moderate  Storms 

180.  A  verified  model  enables  examination  of  the  shore  protection 
functioning  of  alternative  fill  cross  sections  or  fill  templates  at  a  beach 
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nourishment  project  impacted  by  storms .  Such  calculations  must  be  made  to 
perform  cost-benefit  analysis  and  to  estimate  the  minimum  fill  volume  needed 
to  provide  a  specified  level  of  storm  protection  against  erosion  and  inunda¬ 
tion.  In  an  overall  evaluation,  longshore  variability  of  the  protective  beach 
complex  as  well  as  storms  must  be  considered  (van  de  Graaff  1983,  Savage  and 
Birkemeier  1987) . 

181.  Two  approaches  can  be  taken  to  estimate  storm  impact  on  a  beach- 
fill  complex:  the  design- storm  approach  and  the  storm-ensemble  approach.  The 
design  storm  is  either  a  hypothetical  or  a  historical  event  that  produces  a 
specified  storm  surge  hydrograph  and  wave  condition  at  the  project  site. 

Surge  is  a  water  level  rise  caused  by  wind  stress  and  atmospheric  pressure 
variation;  waves  also  produce  a  rise  in  mean  water  level  through  setup.  The 
time  average  (on  the  order  of  an  hour)  of  surge,  wave  setup,  and  tide  is 
called  the  stage.  In  a  stage -frequency  analysis,  the  design  storm  may  have  a 
certain  frequency  of  occurrence  or  return  period  (for  example,  a  100-year 
storm) . 

182.  The  design-storm  approach  is  problematic  for  use  in  beach  erosion 
modeling  because  profile  change  is  sensitive  to  storm  duration,  surge  shape, 
and  wave  height  and  period,  in  addition  to  peak  stage.  The  maximum  water 
level  associated  with  the  surge  of  a  design  storm  may  produce  less  erosion 
than  a  storm  of  lower  surge  but  longer  duration,  or  a  storm  of  lower  surge  but 
higher  waves.  However,  in  the  absence  of  an  extensive  storm  data  base,  use  of 
several  large  historical  storms  (large  in  the  sense  of  having  caused  signifi¬ 
cant  erosion)  to  drive  the  model  is  considered  acceptable.  In  this  approach, 
it  is  recommended  that  the  phasing  of  the  peak  tide  be  varied  with  the  peak 
surge  because  model  results  are  sensitive  to  water  level.  For  example, 
arrival  of  the  peak  surge  during  low  tide  on  a  mesotidal  or  macrotidal  coast 
may  produce  little  or  no  erosion,  whereas  arrival  of  the  peak  surge  at  high 
tide  could  bring  disaster. 

183.  The  most  satisfactory  solution  to  the  problem  of  the  many-to-one 
relation  between  beach  erosion  and  stage  frequency  is  to  use  the  storm- 
ensemble  approach,  i.e.,  calculating  erosion  for  a  large  number  of  storms  and 
keying  the  erosion  to  the  frequency  of  storm  occurrence  (Kriebel  1982; 
Scheffner  1988,  1989).  This  method  yields  an  erosion  or  recession  frequency 
of  occurrence  curve.  In  this  method,  tropical  and  extratropical  storms  must 
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be  treated  independently  because  they  have  different  physical  characteristics. 
Tropical  storms  are  infrequent  events  of  short  duration  and  high  intensity, 
whereas  extratropical  storms  are  more  frequent  and  usually  of  longer  duration 
and  lower  intensity.  The  storm -ensemble  approach  is  recommended  for  project 
design,  although  it  requires  a  storm  data  base  and  is  much  more  computation¬ 
ally  intensive  than  the  design- storm  approach. 

184.  In  this  section,  the  response  of  nourished  beaches  to  a  represen¬ 
tative  tropical  storm  and  an  extratropical  storm  is  calculated  to  examine 
predictions  of  the  model  to  storms  of  differing  waves  and  water  levels  (surge 
and  duration).  The  two  storms  were  synthesized  to  produce  erosion  resulting 
from  a  2-  to  5-year  storm  surge  for  the  mid-Atlantic  Ocean  coast  of  the  United 
States.  The  amount  of  eroded  volume  for  such  events  is  on  the  order  of  20  to 
30  ia?/m  (Savage  and  Birkemeier  1987). 

Design  storms 

185.  Surge  hydrographs  developed  for  the  hypothetical  storms  are  shown 
in  Figure  26.  The  hurricane  surge  has  a  duration  of  approximately  12  hr,  with 
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a  peak  surge  of  2  m,  and  a  duration  above  half  the  peak  surge  (1  m)  of  6  hr. 
The  shape  of  the  hydrograph  follows  an  inverse  hyperbo'ic  cosine  squared.  The 
surge  of  the  northeaster  has  a  duration  of  36  hr,  with  a  peak  surge  of  1  m, 
and  a  duration  above  half  the  peak  surge  (0.5  m)  of  18  hr.  A  cosine -squared 
function  was  used  to  generate  the  shape  of  the  northeaster  surge.  The  peak 
surge  of  the  hurricane  is  higher  because  wind  speeds  in  hurricanes  are  on 
average  greater  than  those  in  northeasters. 

186.  The  time-histories  of  wave  height  and  period  assigned  to  the 
storms  are  shown  in  Figure  27.  Both  have  peak  deepwater  wave  heights  of  5  m, 
which  occur  during  the  time  when  the  respective  surges  are  greater  than  half 
the  maximum.  The  duration  of  high  waves  for  the  northeaster  is  thus  three 
times  that  of  the  hurricane.  Since  the  radius  of  a  northeaster  is  typically 
several  times  greater  than  that  of  a  hurricane,  the  fetch  is  longer,  resulting 
in  longer  wave  periods.  Wave  height  and  period  of  1  m  and  7  sec  were  applied 
for  6.5  days  before  the  start  of  the  storms  to  mold  the  profiles  into  a 
realistic  shape  under  ordinary  waves.  As  the  storm  passed,  wave  height  and 
period  were  changed  to  0.5  m  and  10  sec  to  simulate  long-period  recovery  swell 
wave  conditions.  A  sinusoidal  tide  was  applied  with  a  12-hr  period  and  0.5-m 
amplitude;  a  peak  in  the  tide  occurred  during  the  peak  surge  of  each  storm. 
Beach- fill  shapes 

187.  Following  the  procedure  of  Kraus  and  Larson  (1989),  two  different 
beach-fill  cross  sections  or  templates  were  designed  for  exposure  to  storms. 
The  artificial  berm  had  most  of  the  fill  placed  on  the  beach  and  above  MSL. 

The  as-constructed  berm  extended  horizontally  for  a  distance  of  16  m  at  an 
elevation  of  3  m  and  then  tapered  with  a  1:20  slope  to  join  the  original  beach 
profile  at  a  depth  of  1.4  m  (Figure  28).  This  is  a  common  beach-fill  design 
in  the  United  States  (e.g.,  Jarrett  and  Hemsley  1988).  The  other  fill 
template  is  termed  profile  nourishment,  a  concept  advocated  by  Bruun  (1988). 

In  the  example,  the  material  was  placed  over  the  profile  from  an  elevation  of 
1  to  -2  m  (Figure  29)  in  an  approximation  of  the  existing  profile  shape.  The 
fill  volume  was  the  same  for  each  template  (140  m^/m) . 

188.  Initially,  all  simulations  were  made  assuming  a  0.2-mm  sand  beach, 
with  the  fills  and  the  natural  beach  having  the  same  grain  size.  However, 
simulations  were  also  made  for  fill  grain  sizes  in  0.1-mm  increments  from 

0.2  to  1.0  mm.  In  these  cases,  the  grain  size  was  specified  as  that  of  the 
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Figure  27.  Time-history  of  wave  height  and  period 
for  synthetic  storms 
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fill  over  the  portion  of  the  profile  originally  occupied  by  the  fill  and 
0.2  mm  elsewhere.  Profile  elevations  are  given  with  respect  to  MSL,  and  a 
water  temperature  of  20°  C  was  specified  in  the  model  for  computation  of  the 
sand  fall  speed. 

Profile  change 

189.  Figures  28  through  31  illustrate  predicted  impacts  of  the  two 
storms  on  the  two  nourishment  schemes.  The  curve  labeled  "profile  without 
fill"  gives  a  hypothetical  dune,  beach,  and  subaqueous  equilibrium  profile  for 
reference.  The  construction  profile  configuration  is  denoted  as  "profile  with 
fill,"  whereas  the  "prestorm  profile"  shows  the  profile  after  6.5  days  of 
ordinary  waves.  The  "poststorm  profile"  represents  the  shape  prior  to  the 
recovery  wave  period  and  the  "poststorm  recovery"  profile  displays  the  profile 
after  experiencing  approximately  2  weeks  of  recovery  waves.  The  following 
paragraphs  discuss  prestorm  and  poststorm  profile  development. 

190.  Prestorm.  Prestorm  profiles  ("design  profile"  of  Figure  25)  of 
the  artificial  berm  and  nourished  profile  differ  significantly  in  the  inner 
surf  zone.  A  pronounced  scarp  is  produced  with  the  artificial  berm,  whereas 
the  nourished  profile  experiences  more  subtle  changes  because  it  was  placed  in 
a  near-equilibrium  configuration.  For  both  cross  sections,  a  small  breakpoint 
bar  formed  at  about  the  210-  to  220-m  mark  (measured  from  an  arbitrary 
baseline) .  Sand  was  removed  from  the  foreshore  and  distributed  along  the 
profile  just  beyond  the  2-m-depth  contour.  Thus,  regardless  of  the  initial 
fill  configuration,  the  model  predicts  that  fill  material  will  be  transported 
a  considerable  distance  offshore  to  the  point  of  incipient  wave  breaking  in 
the  process  of  molding  the  surf  zone  profile  to  an  equilibrium  shape. 

191.  Poststorm.  Subaqueous  portions  of  the  poststorm  profiles  are 
quite  similar  in  shape,  being  reworked  by  storm  waves  to  the  same  equilibrium 
shape  independent  of  initial  profile  configuration  and  type  of  storm.  The 
shoreline  position  (0-depth  contour)  actually  moved  seaward  of  its  prestorm 
location,  with  the  accreted  material  supplied  from  the  ordinarily  subaerial 
portion  of  the  profile  that  was  inundated  due  to  storm  surge.  A  small  bar 
formed  at  approximately  the  5-m-depth  contour  under  the  high  storm  waves,  but 
it  is  not  shown  here  to  better  display  changes  near  the  beach. 

192.  An  important  outcome  of  the  predictions  is  that,  under  the  action 
of  simulated  hurricane  and  northeaster  events,  resultant  profile  change  was 


89 


Elevation  (m) 


Elevation  (m) 


Figure  31.  Response  of  nourished  profile  to  northeaster  impact 
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very  similar.  The  use  of  one  typical  storm  descriptor  (e.g.  maximum  stage)  to 
estimate  shoreline  recession  or  eroded  volume  would  have  predicted  greater 
changes  associated  with  the  hurricane,  thus  producing  misleading  results. 

193.  Poststorm  recovery.  In  all  cases,  substantial  accretion  occurred 
near  MSL  during  the  poststorm  recovery  phase.  The  upper  foreshore  of  the 
nourished  profile  experienced  more  accretion  than  the  artificial  berm  cases. 
These  results  are  consistent  with  the  concept  that  a  beach  profile  with  a 
natural  shape  can  best  respond  to  changes  in  incident  waves. 

Eroded  volume  and  contour  change 

194.  Figure  32  shows  the  evolution  of  eroded  volume  above  the  0-  and 
1-m  contours  (and  0.5-m  contour  for  the  profile  nourishment  case).  Eroded 
volume  above  the  0-depth  contour  increases  rapidly  at  the  beginning  of  pre¬ 
storm  ("ordinary")  wave  action,  describing  the  behavior  of  fill  material 
during  initial  profile  adjustment.  In  contrast,  eroded  volume  above  the  1-m 
contour  shows  a  much  slower  increase.  The  profile  nourishment  case  illu¬ 
strates  greater  initial  erosion  during  the  early  stage  of  wave  action,  but 
also  greater  recovery  during  the  poststorm  period.  It  should  be  emphasized 
that  longshore  processes  are  omitted  from  the  present  discussion;  however, 
profile  nourishment  places  a  greater  amount  of  material  into  the  active 
littoral  zone,  thereby  increasing  its  potential  for  transport  alongshore  and 
out  of  the  project  reach. 

195.  The  volume  of  eroded  material  above  the  0-depth  contour  does  not 
show  a  significant  increase  during  the  storms,  changing  from  about  22  to 

27  m^/m  for  erosion  of  the  artificial  berm  during  the  northeaster.  Eroded 
volume  above  the  1-m  contour  abruptly  increases  at  the  start  of  the  storms 
from  about  3  to  17  m^/m  for  the  artificial  berm  impacted  by  the  northeaster. 
The  volume  of  erosion  above  the  0-depth  contour  is  relatively  unchanged 
because  these  moderate  storms  primarily  remove  material  from  the  upper  portion 
of  the  profile  and  redistribute  it  over  the  shallow  nearshore  region. 

196.  The  hurricane  and  the  northeaster  produce  nearly  the  same  amount 
of  erosion,  25  to  30  ra^/m  above  the  0-depth  contour  and  13  to  16  w?/m  above 
the  1-m  contour.  Therefore,  the  longer  surge  duration  of  the  northeaster  was 
approximately  equivalent  in  erosion  capacity  to  the  higher  surge  of  the 
shorter  duration  hurricane.  Eroded  volumes  above  the  0-depth  contour  are 
comparable  with  those  associated  with  2-  to  5-year  return  period  storms 
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Figure  32.  Eroded  volume  above  specific  contours 
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impacting  the  mid-Atlantic  coast  (Savage  and  Birkemeier  1987).  Evolution  of 
the  eroded  volumes  above  MSL  shows  an  approximate  exponential  approach  to  an 
equilibrium  value.  These  trends  are  in  general  agreement  with  those  obtained 
by  Kriebel  and  Dean  (1985),  who  numerically  examined  eroded  volume  and  berm 
recession  as  a  function  of  wave  height,  surge  level,  and  other  parameters. 

197.  It  is  interesting  that  there  is  a  tendency  for  the  hurricane- 
impacted  beach  to  recover  more  rapidly  than  the  northeaster- impacted  beach 
(Figure  32) .  This  is  in  agreement  with  qualitative  observations  of  poststorm 
recovery  along  the  Florida  coast.  An  explanation  lies  in  the  nature  of  the 
storm  hydrograph;  material  is  not  removed  as  far  offshore  during  hurricanes  as 
in  the  case  of  the  longer  duration  and  lower  surge  extratropical  storms.  This 
material  is  deposited  closer  to  shore  and  can  be  moved  back  on  the  beach  more 
rapidly  under  recovery  wave  action. 

198.  Figure  33  displays  the  evolution  of  the  position  of  the  0  and  1-m 
contours.  A  decrease  in  contour  position  indicates  recession  of  the  beach  at 
that  contour.  The  0-depth  contour  for  both  fill  types  and  both  storm  types 
shows  recession  durirg  prestorm  and  storm  periods,  but  begins  to  advance 
before  the  end  of  the  storms  and  prior  to  arrival  of  recovery  waves,  as  the 
surge  subsides  and  the  wave  height  decreases.  The  1-m  contour  shows  no 
recovery  for  the  artificial  berm  because  the  foreshore  is  relatively  steep, 
whereas  the  0.5-m  contour  for  the  nourished  profile  does  show  some  recovery 
since  the  gentler  slope  allows  poststorra  wave  runup  to  build  a  berm.  In 
SBEACH,  berm  formation  and  growth  are  largely  controlled  by  the  elevation 
reached  by  wave  runup  (Larson  1988,  Larson  and  Kraus  1989a),  and  accretionary 
processes  are  diminished  on  steeper  slopes. 

199.  Figure  34  illustrates  eroded  volume  at  the  end  of  the  storm  (prior 
to  recovery  wave  action)  as  a  function  of  beach-fill  grain  size,  As  previ¬ 
ously  mentioned,  median  grain  diameter  of  the  native  beach  was  set  at  0.2  mm, 
and  the  area  in  which  the  fill  was  placed  was  assigned  the  grain  size  of  the 
fill.  This  procedure  does  not  allow  tracking  of  movement  of  the  different 
grain  sizes.  However,  because  surf  zone  sediment  is  usually  sorted  with 
coarser  material  located  higher  on  the  active  profile,  this  simple  procedure 
is  considered  to  provide  a  reasonable  first  approximation  of  the  response  of  a 
natural  beach  of  varying  grain  size. 
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Figure  33.  Movement  of  selected  contours 


Figure  34.  Effect  of  grain  size  on  eroded  volume 
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200.  For  the  particular  synthetic  storms  used,  Figure  34  shows  a  rela¬ 
tively  steep  decrease  in  eroded  volume  as  grain  size  increases  from  0.2  to 
0.4  mm,  with  a  gentle  decrease  thereafter.  This  behavior  results  from  the 
property  of  the  empirically  determined  functional  dependence  of  wave  energy 
dissipation  needed  to  generate  an  equilibrium  profile  of  a  given  grain  size. 
The  property  is  manifested  by  computed  dissipation  rates  that  rise  steeply  in 
the  range  of  0.1  to  0.4  mm  and  then  increase  at  a  lower  rate  with  increasing 
grain  size  (Moore  1982) .  Since  the  rate  of  decrease  in  erosion  is  small 
beyond  0.4  mm  and  the  cost  of  beach  fill  typically  increases  substantially  for 
larger  size  material,  calculations  such  as  those  illustrated  in  Figure  34 
allow  an  evaluation  to  be  made  of  initial  nourishment  volume  and  subsequent 
fill  maintenance  costs.  In  project  planning,  curves  similar  to  those  in 
Figure  34  should  be  developed  with  consideration  of  the  design  profile  and  the 
storm  characteristics  of  the  area,  using  the  storm- ensemble  approach. 

Beach- Fill  Response  with  a  Seawall 

201.  Another  aspect  of  the  simulation  was  to  display  the  capabilities 
of  the  model  to  predict  profile  evolution  of  a  fill  fronting  a  seawall.  A 
schematic  storm  event  was  developed  similar  to  the  previous  two  synthetic 
storms  and  consisted  of  a  period  of  low  waves  when  profile  adjustment  took 
place,  a  3 -day  storm  event,  and  a  period  of  recovery  waves.  The  total 
simulation  time  was  30  days,  and  comparisons  were  made  for  a  beach  with  and 
without  fill.  The  two  fill  types  studied  were  again  an  artificial  berm  and 
profile  nourishment,  with  cases  examined  using  median  grain  diameters  of 

0.25  and  0.40  mm.  Although  wave  reflection  from  the  seawall  is  not  considered 
in  the  simulations,  it  is  recognized  that  this  process  alters  surf  zone  hydro¬ 
dynamics  and  may  affect  profile  development. 

202.  The  first  21  days  of  the  simulation  consisted  of  constant  waves 
and  water  level,  with  a  wave  height  of  0.5  m  and  a  wave  period  of  8  sec. 

These  values  were  assumed  to  represent  typical  wave  conditions  on  an  open- 
ocean  Atlantic  coast  and  were  used  to  produce  a  realistic  equilibrium  profile. 
Under  these  conditions,  beach  change  can  be  mildly  erosional  or  accretionary, 
depending  on  the  grain  size  and  initial  profile  shape.  At  day  21  (" Prestorm" 
in  Figure  35)  a  storm  was  imposed  that  lasted  for  3  days.  The  average  surge 
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Figure  35.  Time -history  of  waves  and  water  level 


water  level  during  the  storm  followed  a  Gaussian  form  and  reached  a  maximum 
elevation  of  3  m  above  MSL.  The  average  wave  height  was  given  a  sinusoidal 
variation  during  the  storm,  having  a  peak  of  2.5  m  at  the  peak  of  the  surge 
("Midstorm").  Wave  period  during  the  3-day  event  varied  between  8  and  12  sec, 
with  12 -sec  waves  occurring  at  the  time  of  the  maximum  wave  height.  At  the 
end  of  the  storm  ("Poststorm"),  the  surge  was  considered  nonexistent,  and 
waves  arrived  as  swell  with  a  height  of  0.5  m  and  a  period  of  16  sec. 

203.  In  the  following  sections,  adjustment  of  the  profile  to  the  30-day 
design  sea  conditions  is  examined  with  SBEACH  for  three  situations:  an 
original  (existing)  condition,  a  fill  of  artificial  berm- type,  and  a  fill 
carried  out  using  the  profile  nourishment  technique.  The  simulations  were 
done  using  parameters  from  the  calibrated  model  with  a  time -step  of  20  min  and 
a  space  step  of  5  m. 

Existing  conditions 

204.  The  profile  configuration  representing  the  beach  prior  to  nourish¬ 
ment  activity  is  shown  in  Figure  36.  A  2-m-high  foreshore  is  backed  by  an 
infinitely  high  seawall,  although  its  crown  is  truncated  at  a  2-m  elevation  in 
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Original  Beach  (0.25-mm  Sand) 


this  and  subsequent  figures.  The  original  profile  decreases  linearly  to  the 
0.5-m-depth  contour  with  two  different  slopes;  the  break  in  slope  is  located 
at  an  elevation  of  1  a.  Seaward  of  the  segmented  linear  profile,  an  equili¬ 
brium  shape  is  used  into  deeper  water. 

205.  Simulations  showed  that  the  0.25 -mm- sand  beach  eroded  on  the 
foreshore  during  the  first  21  days  of  low  waves,  and  a  small  bar  formed  with  a 
crest  located  about  120  m  from  the  baseline.  In  contrast,  the  0.40-mm-sand 
beach  accreted  slightly  on  the  foreshore,  indicating  that  the  transport  was 
directed  onshore.  At  midstorm,  two  bars  formed  offshore,  the  berm  was 
completely  submerged,  and  the  beach  eroded  considerably  at  the  seawall  for 
both  sand  sizes.  The  inner  bar  was  created  at  the  beginning  of  the  surge, 
whereas  formation  of  the  outer  bar  corresponded  to  the  peak  surge  and  maximum 
wave  height.  Poststorm  profile  shape  showed  that  the  berm  had  eroded  further 
at  the  seawall,  slightly  more  so  for  the  0.25 -mm- sand  beach  than  for  the 
0.40-mm  sand  beach,  and  a  single  large  bar  appeared  on  both  profiles.  The  bar 
on  the  finer  sand  beach  is  broader  and  lower  in  elevation  with  respect  to  the 
initial  profile  than  the  coarser  sand  bar.  It  is  interesting  to  note  that  the 
position  of  the  poststorm  shoreline  is  at  about  the  same  location  (approxi¬ 
mately  75  m  from  baseline)  for  both  beaches;  however,  sand  moved  farther 
offshore  on  the  finer  sand  beach,  and  its  subaerial  section  is  more  deflated. 

206.  Lower  steepness  waves  arriving  during  the  7-day  recovery  period 
passed  over  the  storm  bars  on  both  beaches,  causing  little  change  in  bnr 
shape.  During  the  recovery  phase,  swell  broke  in  shallower  water,  creating  a 
small  secondary  bar  inshore  on  both  profiles.  The  inner  bar  was  formed  by 
onshore  transport,  producing  a  deep  trough  between  the  storm  bar  and  inshore 
bar.  Both  beaches  accreted  to  the  limit  of  runup. 

Artificial  berm  fill 

207.  Simulation  results  for  the  artificial  3-m-high  berm  are  given  in 
Figure  37.  A  beach  fill  of  volume  85  rn^/rn  with  the  same  grain  size  as  the 
existing  beach  was  placed  mainly  on  the  subaerial  portion  of  the  profile.  The 
fill  extended  horizontally  from  the  seawall  for  20  m,  tapered  down  for  several 
meters,  and  then  extended  offshore  with  slope  1:10  to  meet  the  equilibrium 
profile  at  a  depth  of  approximately  1  m. 

208.  Twenty-one  days  of  prestor  ;aves  eroded  the  toe  of  the  0.25-mra 
artificial  sand  berm  and  created  a  small  bar,  whereas  the  0.40-mm-sand  be-*ch 
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Figure  37.  Artificial  berm  simulations 
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configuration  was  effectively  stable  under  the  mild  waves.  Berm  elevation  for 
both  sand  sizes  was  lowered  during  the  storm,  but  not  to  the  extent  of  the 
original  beach  case  (Figure  36) .  The  artificial  berm  provided  substantial 
protection,  but  should  have  been  higher  to  prevent  waves  from  reaching  the 
seawall  for  this  extreme  storm  surge.  The  prestorm  bar  was  located  closer  to 
shore  than  in  the  existing  case  simulation.  Poststorm  bars  for  both  sand 
sizes  were  located  in  essentially  the  same  position  as  for  the  corresponding 
existing  beach  cases. 

209.  The  recovery  cross  sections  for  both  sand  sizes  were  considerably 
different  from  the  original  beach  cases  because  of  the  greater  amount  of  sand 
in  the  system.  Although  significant  recovery  took  place  for  the  0.25-mm-sand 
beach,  sand  was  effectively  trapped  in  the  storm  bar  and  could  not  contribute 
to  poststorm  beach  recovery.  Eroded  material  for  the  coarser  sand  beach  was 
deposited  closer  to  shore  and,  during  the  recovery  phase,  the  entire  storm 
profile  translated  shoreward.  Although  a  limited  amount  of  foreshore  recovery 
took  place  with  the  0.40-mm  sand,  a  significant  volume  of  material  remained 
near  to  shore. 

Profile  nourishment 

210.  The  profile  nourishment  case  is  depicted  in  Figure  38,  where  the 
major  amount  of  material  is  placed  along  the  nearshore  portion  of  the  profile 
in  accordance  with  the  concept  that  a  beach  can  best  protect  itself  if  in  a 
natural  shape.  The  same  volume  of  fill  as  in  the  artificial  berm  case 

(85  m^/m)  was  placed  over  the  profile  in  an  equilibrium  shape  and  tapered  at 
the  seaward  end.  This  fill  template  resulted  in  positioning  of  the  shoreline 
approximately  10  m  seaward  of  the  artificial  berm  case. 

211.  Simulation  results  show  that  subaerial  erosion  was  reduced 
significantly  as  compared  with  existing  conditions,  even  though  the  entire 
profile  was  submerged  for  much  of  the  storm.  Two  bars  developed  again,  but 
the  inner  bar  ("Midstorm")  for  the  0.25 -mm  sand  generated  during  the  rising 
surge  was  located  shoreward  of  the  prestorm  bar.  Bar  development  (with 
respect  to  the  initial  profile)  for  this  grain  size  was  more  pronounced  than 
that  for  the  artificial  berm  case,  whereas  profile  evolution  with  0.40-mm  sand 
showed  less  bar  development.  The  profile  nourishment  template  showed  greater 
retention  of  sand  above  MSL  than  the  artificial  berm  case,  and  less  material 
was  redistributed  along  the  profile.  However,  because  of  its  low  elevation, 
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Figure  38.  Profile  nourishment  simulations 


the  profile  nourishment  template  exposed  the  wall  to  more  wave  action  during 
the  surge  than  did  the  artificial  berm  of  the  same  volume. 

Discussion  of  Beach- Fill  Simulations 

212.  The  numerical  model  allowed  realistic  time -dependent  calculations 
to  be  made  of  beach- fill  change  in  terms  of  readily  available  engineering 
data.  Model  performance  was  realistic  regarding  storm- induced  beach  erosion, 
showing  correct  dependencies  on  surge  and  wave  input,  as  well  as  initial 
profile  shape.  Storm- induced  beach  and  dune  erosion  cannot  be  uniquely 
specified  through  a  single  storm-related  parameter  such  as  maximum  stage. 

This  result  demonstrates  the  limited  usefulness  of  the  design  storm  approach. 
Qualitative  features  of  the  posts torm  recovery  phase  were  also  reproduced,  but 
further  work  is  required  to  accurately  simulate  recovery. 

213.  The  model  can  be  used  to  judge  the  relative  behavior  of  various 
beach- fill  cross  sections  exposed  to  ordinary  and  extreme  waves  for  time 
intervals  on  the  order  of  days  to  weeks.  Fill  placed  on  the  upper  beach  was 
calculated  to  move  offshore  to  the  region  of  incipient  wave  breaking,  in 
general  agreement  with  the  initial  movement  of  beach-fill  material  (Hansen,  in 
preparation).  Also,  a  limited  number  of  calculations  indicated  that  it  may 
not  be  cost-effective  to  use  beach  fill  with  a  median  grain  diameter  much 
greater  than  0.4  mm  for  storm  protection,  owing  to  the  typically  greatly 
increased  cost  and  the  declining  benefit  in  decreased  volume  of  eroded 
material.  Design  alternatives  for  specific  situations  (fill  grain  size,  fill 
cross  section,  and  storm  climate)  could  be  evaluated  with  the  model. 

214.  The  model  is  based  on  relationships  for  cross -shore  sand  transport 
produced  by  breaking  waves,  implying  that  longshore  sand  transport  is  con¬ 
sidered  negligible.  Thus,  caution  must  be  exercised  in  applying  the  model  and 
interpreting  results  if  gradients  in  longshore  transport  exist.  Longshore 
transport  is  expected  to  be  a  major  contributing  process  to  beach  change  over 
long  periods,  for  example,  seasons  to  a  year.  For  such  time  frames,  the 
present  model  should  be  combined  with  other  predictive  methodologies,  such  as 
a  shoreline  change  model  (e.g.  Kraus  1983;  Kraus,  Hanson,  and  Harikai  1984; 
Hanson  1989;  Kraus  1989),  which  compute  the  time  rate  of  change  of  shoreline 
position  based  on  estimates  of  wave -induced  longshore  sand  transport. 


103 


PART  VII:  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 


215.  The  numerical  model  of  beach  profile  change  SBEACH  was  developed 
using  large  wave  tank  data  and  was  verified  with  field  data,  as  reported  in 
previous  work  (Larson  1988,  Larson  and  Kraus  1989a).  In  the  present  report, 
additional  model  verification  was  made  using  field  data  on  berm  and  dune 
erosion  from  the  US  east  and  west  coasts.  Overall  model  performance  for 
simulating  the  erosional  phase  of  a  storm  event  was  satisfactory,  whereas  the 
recovery  phase  following  a  storm  event  was  less  well  described  and  essentially 
only  in  qualitative  agreement  with  the  measured  profile  evolution. 

216.  A  high-quality  data  set  describing  storm- induced  erosion  off  the 
coate  of  New  Jersey  was  used  to  investigate  the  capability  of  the  model  to 
simulate  berm  and  dune  erosion.  SBEACH  predicted  eroded  volumes  and  contour 
retreat  in  agreement  with  what  could  be  confidently  inferred  to  occur  during 
the  erosional  part  of  the  storm  event.  When  the  storm  waned  and  waves  became 
accretionary,  a  large  sand  ridge  formed  on  the  foreshore  of  a  size  that  was 
not  predicted  by  the  model.  However,  qualitative  features  of  the  recovery 
were  reproduced  by  the  model,  including  development  of  profiles  with  different 
shapes  according  to  initial  conditions. 

217.  The  capability  of  the  model  to  describe  profile  change  on  the  US 
west  coast  was  evaluated  with  data  available  from  minor  storms  at  Torrey  Pines 
Beach,  California.  Contour  retreat  on  the  foreshore  was  well  simulated  along 
the  general  profile  shape.  The  amount  of  erosion  predicted  by  the  model  was 
also  in  good  agreement  with  field  measurements.  Similar  quantities  of  calcu¬ 
lated  and  measured  erosion  were  obtained  by  increasing  the  water  level  by 

0.3  m  over  that  measured  at  the  gage  located  offshore  in  9.1  m  of  water.  This 
increase  can  be  ascribed  to  a  number  of  factors  such  as  locally  generated  wind 
setup  landward  of  the  gage,  wave  setup,  and  differences  in  atmospheric 
pressure.  Further  theoretical,  laboratory,  and  field  research  is  needed  to 
better  quantify  hydrodynamic  and  sand  transport  processes  in  the  swash  zone, 
as  model  predictions  are  sensitive  to  changes  in  water  level  nearsnore. 

218.  SBEACH  was  applied  to  simulate  the  short-term  fate  of  beach  fill  by 
examining  initial  adjustment  of  the  fill  and  its  response  to  storms.  In  these 
hypothetical  cases,  beach-fill  cross  section  and  grain  size  were  varied,  and 
the  response  of  the  fill  to  synthetic  storm  events  was  investigated.  Synthetic 


storms  corresponding  to  tropical  (hurricane)  and  extratropical  (northeasters) 
events  were  chosen  as  being  typical  for  the  US  east  coast.  Two  fill  cross 
sections  were  evaluated,  the  artificial  berm  and  profile  nourishment. 

219.  Simulations  of  fill  response  to  storm  events  showed  the  importance 
of  peak  surge  level,  surge  duration,  and  wave  conditions.  For  example,  a  high 
surge  with  a  low  duration  produced  as  much  erosion  as  a  low- surge  storm  with  a 
longer  duration.  The  profile  nourishment  fill  design  showed  less  change  for 
the  design  storm  than  the  artificial  berm  of  the  same  volume,  suggesting  that 
the  nourished  profile  was  more  in  accordance  with  the  natural  equilibrium 
profile  and  minimized  profile  change  as  the  wave  and  water  level  conditions 
changed.  This  assumption  was  further  substantiated  by  the  stronger  recovery 
shown  by  the  profile -nourished  beach.  Model  simulation  also  indicated  that 
little  gain  in  fill  stability  (erosion  protection)  was  achieved  during  a  storm 
if  the  fill  grain  size  exceeded  0.40  mm.  Overall,  the  numerical  model 
provided  reasonable  results  for  simulating  the  response  of  beach  fills  and 
could  be  used  in  a  cost-benefit  analysis  for  determining  the  optimal  fill 
design  and  protection  level.  However,  it  is  recommended  that  the  model  be 
verified  prior  to  use  in  planning  and  design  with  data  from  the  project  site 
or  a  neighboring  representative  site. 

220.  Conditions  under  which  the  numerical  model  is  expected  to  give 
reliable  results  may  be  summarized  as  follows: 

4.  Variations  in  water  level  and  breaking  waves  are  the  main 
causes  of  sand  transport  and  thus  profile  change. 

Longshore  effects  are  negligible,  that  is,  the  incident  waves 
approach  normal  to  the  shoreline,  or  differentials  in  longshore 
transport  are  minor  (influence  of  structures  blocking  longshore 
transport  is  small,  and  the  shoreline  is  straight). 

£.  Median  grain  diameter  on  the  profile  is  in  the  interval  0.1  to 
2.0  mm  and  reasonably  uniform  across -shore. 

221.  Experience  in  simulating  field  situations  has  identified  several 
areas  requiring  basic  research  and  model  refinement.  The  most  important  area 
for  further  research  is  the  description  of  onshore  sand  transport  and  berm 
development.  The  poststorm  recovery  phase  must  be  reproduced  to  properly 
estimate  storm  impacts,  since  field  surveys  may  be  made  after  some  recovery 
has  already  taken  place. 
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222.  SBEACH  and  similar  beach  erosion  models  are  sensitive  to  changes 
in  water  level.  Considerable  effort  is  required  to  make  progress  in  predict¬ 
ing  water  levels  at  the  coast  in  an  economical  and  complete  fashion.  Pheno¬ 
mena  to  be  studied  and  incorporated  into  a  coastal  engineering  methodology 
compatible  with  profile  change  modeling  include  wind  setup,  infragravity  wave 
motion,  and  contributions  from  differentials  in  atmospheric  pressure. 

223.  In  the  present  version  of  SBEACH,  the  same  transport  formula  is 
used  independent  of  transport  direction.  This  description  is  probably  accep¬ 
table  in  the  initial  stage  of  beach  recovery  when  waves  break  offshore  and  a 
well-established  surf  zone  exists.  However,  as  material  is  transported 
onshore  and  waves  break  closer  to  shore,  the  width  of  the  surf  zone  is 
reduced,  and  the  character  of  wave  breaking  might  change  abruptly  from 
plunging  to  surging.  This  development  during  the  recovery  phase  is  not 
described  satisfactor 'ly  by  the  model.  Also,  onshore  transport  in  deeper 
water  outside  the  breax  point,  which  takes  place  over  long  periods,  lacks  a 
good  description  in  the  model.  Finally,  regarding  the  transport  rate  formula¬ 
tion,  only  two  types  of  cross -shore  distributions  are  presently  allowed, 
purely  onshore  or  purely  offshore  transport.  During  times  of  changing  wave 
conditions  and  mixed  waves,  more  complex  transport  rate  distributions  occur, 
particularly  if  the  profile  is  nearly  in  equilibrium  with  the  waves. 

224.  It  is  possible  to  specify  different  grain  sizes  (equilibrium 
energy  dissipation)  according  to  location  along  the  profile  in  the  model,  but 
since  no  tracking  of  individual  grains  is  performed,  the  description  is  rather 
crude.  Thus,  it  is  only  realistic  in  the  present  model  to  simulate  beach 
conditions  involving  a  few  distinct  grain  sizes.  By  keeping  track  of  the 
grain  movement  and  assigning  a  probability  distribution  to  the  grains  in  each 
calculation  cell,  it  would  be  possible  to  improve  this  aspect  of  the  model. 

225.  Future  research  for  refining  and  extending  the  capabilities  of 
SBEACH  will  be  directed  toward  improved  understanding  and  modeling  of  sand 
transport  associated  with  wave  runup,  dune  overtopping,  and  accretionary 
processes  on  the  foreshore;  study  of  the  representation  and  action  of 
irregular  waves;  and  collection  of  high-quality  field  data.  At  present,  model 
reliability  is  limited  by  necessity  to  calibrate  the  model  at  the  project  site 
or  neighboring  sites.  Field  data  collected  immediately  before,  during,  and 


after  storm  events  causing  significant  change  in  beach  profile  shape  will 
allow  examination  and  improvement  of  individual  model  algorithms. 

226.  Values  of  the  calibration  parameters  determined  for  the  east  and 
west  coast  case  studies  can  be  viewed  as  nominal  values  used  in  model  applica¬ 
tions.  Sensitivity  tests  and  variation  of  input  parameters  must  be  made  to 
obtain  a  range  of  predictions  in  which  the  actual  storm- induced  beach  change 
may  be  expected  to  fall. 
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APPENDIX  A:  NOTATION 


A  Shape  parameter  for  equilibrium  beach  profile,  m1/3 
C  Wave  phase  speed,  m/sec 

Cg  Wave  group  speed,  m/sec 

CQ  Wave  group  speed  in  deep  water,  m/sec 

d  Total  water  depth,  m 

D50  M'  *ian  grain  diameter,  mm 

D  Wave  energy  dissipation  per  unit  water  volume,  N-m/m3-sec 
De q  Equilibrium  wave  energy  dissipation  per  unit  water  volume,  N-m/m3-sec 
E  Wave  energy  density,  N-m/m2 

Es  Stable  wave  energy  density,  N-m/m2 

F  Wave  energy  flux,  N-m/m- sec 

Fg  Stable  wave  energy  flux,  N-m/m-sec 

g  Acceleration  due  to  gravity,  m/sec2 

h  Still -water  depth,  m 

H  Wave  height,  m 

Hjj  Breaking  wave  height,  m 

Hjjq  Energy-based  significant  wave  height,  m 

H0  Deepwater  wave  height,  m 

Hrms  Root-mean- square  wave  height,  m 

Hso  Deepwater  significant  wave  height,  m 

H  Mean  wave  height,  m 

H0  Mean  deepwater  wave  height,  m 

K  Sand  transport  rate  coefficient,  mVN 

L  Wavelength,  m 

L0  Deepwater  wavelength,  m 

M  Dimensionless  empirical  coefficient 

n  Ratio  of  wave  group  speed  and  wave  phase  speed 
q  Net  cross -shore  sand  transport  rate,  m3/m-sec 
Sxx  Radiation  stress  component  directed  onshore,  N/m 
t  Time,  sec 

T  Wave  period,  sec 

w  Sand  fall  speed,  m/sec 

x  Cross -shore  coordinate,  m 


Al 


y  Longshore  coordinate,  m 

ZR  Height  of  active  subaerial  profile,  m 

P  Local  beach  slope  seaward  of  break  point 

y  Ratio  between  wave  height  and  water  depth  at  breaking 
F  Empirical  stable  wave  height  coefficient 

A  Change  in  quantity 

e  Slope  related  sand  transport  rate  coefficient,  m2/sec 
•q  Mean  water  surface  elevation  (wave  setup  or  setdown) ,  m 

k  Wave  decay  coefficient 

A  Spatial  decay  coefficient,  m*1 

0  Wave  angle  with  respect  to  bottom  contours 
p  Density  of  water,  kg/m3 

i|>  Local  profile  slope 

i|>iy  Angle  of  initial  yi-Ad 

i|>r.  Residual  angle  after  shearing 

5  Surf  similarity  parameter 

Subscripts  and  Superscripts: 
b  Breaking  condition 

eq  Equilibrium  condition 

o  Deepwater  condition 

i  Variable  at  grid  point  or  cell  i 

k  Variable  at  time  step  k 

N  Number  of  cells  where  avalanching  occurs 

p  Plunge  point 

r  Runup  limit 

s  Stable  condition  for  depth- limited  wave 

z  End  of  surf  zone 

1  Specific  value  of  variable  in  zone  I 

2  Specific  value  of  variable  in  zone  XI 

'  Denotes  the  smoothed  transport  rate 


A2 


